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ABSTRACT 

We present a Giant Molecular Cloud (GMC) catalog toward M33, containing 71 GMCs in 

total, based on wide field and high sensitivity CO(J = 3 — 2) observations with a spatial resolution 
of 100 pc using the ASTE 10 m telescope. Employing archival optical data, we identify 75 young 
stellar groups (YSGs) from the excess of the surface stellar density, and estimate their ages by 
comparing with stellar evolution models. A spatial comparison among the GMCs, YSGs, and Hll 
regions enable us to classify GMCs into four categories: Type A showing no sign of massive star 
formation (SF), Type B being associated only with Hll regions, Type C with both Hn regions 
and < IGMyr-old YSGs and Typc-D with both Hii regions and 10 30 Myr YSGs. Out of 65 
GMCs (discarding those at the edges of the observed fields), 1 (1%), 13 (20%), 29 (45%), and 22 
(34%) are Types A, B, C, and D, respectively. We interpret these categories as stages in a GMC 
evolutionary sequence. Assuming that the timescale for each evolutionary stage is proportional to 
the number of GMCs, the lifetime of a GMC with a mass > 10^ M© is estimated to be 20-40 Myr. 
In addition, we find that the dense gas fraction as traced by the CO(J = 3 — 2)/C0(J =1 — 0) 
ratio is enhanced around SF regions. This confirms a scenario where dense gas is preferentially 
formed around previously generated stars, and will be the fuel for the next stellar generation. In 
this way, massive SF gradually propagates in a GMC until gas is exhausted. 

Subject headings: Galaxies: individual (M33) — ISM: clouds — ISM:molecules — Submillimeter: ISM 
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1. Introduction 

Observations in our Galaxy and in the Lo- 
cal Group of Galaxies (LGGs) have revealed that 
Giant Molecular Clouds (GMCs) have masses of 
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about 10^ to 10^ M(7) and sizes of 50 to several hun- 



mura 



dred of parsecs (e.g., Blitz 1993 Fukui & Kawa- 



2010 1, which are often found associated with 



Hii regions and OB associations (Massey 20031, 
tracers of massive star formation (SF). Massive 
stars, which are commonly found clustered in OB 



associations (Lada & Lada 2003), are short-lived 



(less than a few lOMyr), and thus spend most of 
their lifetime within the cluster in which they were 
formed. This suggests that GMCs are the princi- 
pal sites of massive SF. However, how SF occurs 
in a CMC, the so-called "GMC evolution" , has re- 
mained poorly understood because it has been dif- 
ficult to perform molecular gas observations with 
a combination of high resolution, sensitivity and 
wide field to enable studies of a sufficiently large 
sample of GMCs. Such studies in LGGs have the 
potential to reveal the density distribution and 
kinematics at GMC scales, as well as its relation 
to ongoing SF. 

Several studies focusing on GMC evolution in 
the Large Magellanic Clouds (LMC) have been 



performed (e.g., Fukui et al. 


1999 


Mizuno et al. 


2001 


|Yamaguchi et al. 


2001 


Kawamura et al. 


2009 


). The GMCs in the LMC are classified into 



three types according to their associated massive 
SF activities: Type i showing no signature of mas- 
sive SF, Type ii being associated with relatively 
small Hii regions, and Type ill with both Hii 
regions and young (less than lOMyr old) stellar 



clusters (Yamaguchi et al. 2001 Kawamura et al. 



2009 1 . While in general the positions of young 
clusters are correlated with GMCs, the clusters 
older than 10 Myr have a weaker or no correlation 
with the GMCs (IKawamura et al.l 120091). This 



classification was interpreted to be a template for 
the GMC evolutionary stages, and the typical life- 
time of GMCs with masses larger than 5 x 10^ Mq 



is estimated to be 20-30 Myr (Blitz et al. 2007 



Kawamura et al. 2009). Analyses of the kinetic 
temperature and density of several GMCs in LMC 
shows that they increase generally as the GMCs 



evolve ( 


Minamidani et al. 2008 


Pineda et al. 2008 


Mizuno et al.|2010 


Minamidani et al.|2011 


1. This 



confirms that the denser and warmer molecular 
gas is directly linked to SF activities. Because 
these studies are only limited to several specific 
GMCs in one irregular galaxy, it is still necessary 
to investigate the properties of GMCs in the dif- 
ferent environments present in other galaxies, in 



order to complete our current knowledge of tem- 
plates for GMC evolution. 

M33 is one of the best laboratories to study 
GMC evolutions under different environments due 



to its proximity (distance of 840 kpc; Freedman 



Freedman 


51°; 


Cor- 



et al. 2001), favorable inclination (i 
belli & Salucci [2000 ), as well as the existence of 



spiral structure and massive SF regions. M33 is 
the second closest spiral galaxy after M31, but 
M33 is less inclined, which is ideal to resolve gas 
components or individual stars with little contam- 
ination along the line of sight. Also, studying 
GMCs over a disk galaxy enable us to compare all 
GMCs at the same distance, unlike in our Galaxy. 

Surveys of ^^CO{J = 1-0) and ^^CO{J = 2-1) 
emission in M33 to investigate the relation be- 
tween molecular gas and massive SF have been 



conducted by several authors (Engargiola et al. 


20031 Heyer et al. 2004| Rosolowsky et al.||2007| 


Gratier et 


al. 2010| Tosaki et al. 2011 Gratier 


et al. 


2012 


). Engargiola et al. 


(|2003D have found 



that more than two-thirds of the GMCs in M33 
are associated with Hii regions within 50 pc, which 



is similar to the proportion in LMC ( jYamaguchi 



et al. 2001 Kawamura et al. 2009). On the 



other hand, the correlation between the molecu- 
lar gas surface densities traced by CO(J = 1 — 
0) emission and that of SF rate (SFR), the so- 
called Kennicutt-Schmidt Law (hereafter K-S law; 
Schmidt||l959 Kennicut"t||1989 ) in a pixel-to-pixel 



analysis is found to be weaker at 80 pc resolutions 



(Onodera et al. 2010) than at lower resolutions 
(200 pc; [Heyer et al.|20"04l ). This suggests that the 
GMCs traced by CO(J =1-0) are weU corre- 
lated with the massive SF sites, but their peaks 
are offset from each other. 

Recently, the K-S law using C0( J = 3 — 2) in 
M33 has shown a tight correlation at 100 pc res- 
olution, unlike CO(J =1 — 0) (Onodera et al. 
2013, PAS J, accepted). In fact, with higher reso- 



lution data, Tosaki et al. (2007) have found that 
the CO(J = 3 — 2) emission is spatially better cor- 
related to massive SF sites than CO (J =1 — 0) in 
the most luminous giant Hii region (GHR) in M33, 
NGC 604. This supports that the dense and warm 
gas traced by C0( J = 3 — 2) is more closely linked 
to the SF sites, as previously argued in four Virgo 
galaxies in Wilson et al. (2009). These results 



are interpreted as higher kinetic temperatures and 
densities in GMCs found close to GHRs, in agree- 
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ment with large velocity gradient (LVG) analysis 
with multi-J CO transition data by [Wilson et al.| 
(|1997D andlMuraoka et al.l (120121). 



As part of the Nobeyama Radio Observatory 
(NRO) M33 All-disk survey of Giant Molecular 



Clouds (MAGiC) project (Tosaki et al. 2011), here 



we present a catalog of M33 GMCs using new 
CO(J = 3 — 2) maps with 25" (corresponding 
to 100 pc) resolution. In order to investigate the 
relation of the GMCs with massive SF, we also 
present a "young stellar group (YSG)" catalog. 
In this work, YSGs can be clusters (typical sizes 
of 15 pc), OB associations (typical sizes of ^15- 
100 pc), and star complexes (sizes of a few hun- 
dred pc) (e.g., |Battineni et al.|1996l|Ivanov|2005 1. 
There are a number of catalogs for clusters, OB as- 
sociations, and star complexes in M33 published 



so far (e.g., Sarajedini & Mancone||2007 hereafter 


SM; IPark & Lee 2007 


San Roman et al.||2010| 


Humphreys & Sandage 


1980 


Ivanov||2005|. How- 



ever, the number of star clusters whose ages have 
been estimated is limited, and the methods to de- 
rive the ages are not homogenous (SM and refer- 
ences therein; Park et al. 2009, hereafter PPL; [San] 



Roman et al.|2009 hereafter SSGH). We also aim 
to estimate their ages in a homogeneous manner 
using the stellar photometry catalog provided by 



Massey et al. (2006) 



The outline of this paper is as follows. The 
observations and data reduction are summarized 
in Section [2] In Section [3] we explain the gen- 
eral procedure for identifying GMCs, as well as 
the identification of YSGs and the estimation of 
their ages. We classify GMCs as a function of the 
age of YSGs associated with them. We also quan- 
tify the SF activities and dense gas fractions in 
the GMCs using the extinct ion- corrected Ha data 
and the CO(J = 3 - 2)/C0(J =1-0) line ra- 
tio, i?3_2/i-0i E^nd then show the variety of phys- 
ical states of the identified GMCs. In Section |4j 
we discuss the relationship between the properties 
of these clouds and their evolutionary stages. Fi- 
nally, we interpret our results in the context of a 
continuous SF in GMCs. The summary of this 
paper is provided in Section [5] 

2. Observations and Data Reduction 

We describe the data that is used throughout 
this paper, including the CO (J = 3 — 2) data to 



identify the GMCs, the CO(J =1-0) data to 
derive i?3_2/i-0: and the optical and infrared data 
to identify the YSG and Hii regions in M33. 

2.1. ASTE CO( J = 3 - 2) Data 

We chose eight regions covering the northern 
and southern spiral arms as well as the galaxy cen- 
ter, where the CO(J =1 — 0) emission is promi- 
nent ([Tosaki et aL][20TI|). The CO(J = 3-2) 



(rest frequency: 345.796 GHz) mapping observa- 
tions were performed between June and Novem- 
ber 2011, using the Atacama Submillimeter Tele- 
scope Experiment (ASTE; jEzawa et al 



2008 



2004 



Kohno 2005) 10-m dish equipped with the 
345-GHz side band separating (2SB) SIS receiver 
CATS345 (CArtridge-Type Sideband-separating 
receiver for ASTE 345-GHz band; |Inoue et al] 
2008[). An XF-Type Digital spectrometer MAC 



(Sorai et al. 2000) was used to cover a velocity 
width of 445 kms~^ with a velocity resolution of 
2.5kms-i at 345 GHz. The On-The-Fly (OTF) 



mapping technique (Mangum et al. 2007), based 



on scanning smoothly and rapidly along one di- 
rection across a rectangular map, was employed to 
obtain the CO(J = 3 — 2) data. The observed area 
covers in total about 80 % of the CO(J = 1-0) 



based molecular gas mass in Tosaki et al. (2011). 

The observations and the data calibration were 
performed for each region as follows. The half 
power beam width (HPBW) of the ASTE 10- 
m telescope is 22" at 345 GHz. The chopper- 
wheel technique was employed to calibrate the 
antenna temperature and the final data has 
been converted into units of main beam bright- 
ness temperature (Tmb = T'^/?7mb), where rjynh is 
the main beam efficiency, measured to be 0.6 ±0.1 
using the standard source Orion KL at (aj2000; 
i5j2ooo)=(05'^35'"14!5, 05°22'29'.'6). Hereafter, all 
the CO intensity measurements are specified in 
Tmb- The typical system temperatures in a single 
side band were usually less than 200 K. The ab- 
solute pointing accuracy was verified by observing 
the C0( J = 3—2) emission of o-Cet at 1-1.5 hr 
intervals, and it was kept better than 4" rms. The 
absolute intensity stability was also monitored us- 
ing o-Cet during the different observing runs and 
was found to be stable within ~ 16%. 

The data reduction was carried out with the 
Nobeyama OTF Software Tools for Analysis and 
Reduction package (NOSTAR, Sawada et al. 
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2008). The "scanning noise" was removed by 



combining scans using the PLAIT algorithm as 



described by Emerson & Graeve ( 1988 ) . The data 



have been convolved with a Bessel-Gauss function 
in the spectral domain to create the data cube. 
The final grid spacing is S'.'O and the angular res- 
olution is 25". The final data cubes for the eight 
regions are characterized by an rms noise of dch = 
16-32 mK in a velocity resolution of 2.5kms~^, 
achieved after 87 hr of total integration. We 
briefly summarize in Table [l] the observational 
and data reduction parameters for the eight ob- 
served regions, including their central positions, 
sizes, integration times, and rms for a 2.5kms~^ 
channel. 

Figure [l] shows the obtained integrated in- 
tensity map over the velocity range Vlsr = 
— 274kms~"'^ to — 94kms~^, where emission is > 



2o'chj overlaid on the Ha image (§ 2.3). The rms 
noise level of the integrated intensity map (CTmom) 
is 0.38 K km s^^ on average. The noise level of 
the CO(J = 3 — 2) integrated intensity map is 
calculated using tTmon = cFdn^NSv, where N is 
the number of integrated channels and 5v is the 
velocity resolution of a channel (2.5 kms^^). The 
observed fields cover the molecular disk up to 
the galactic radius of ~5kpc. Note that we also 
present CO integrated intensity maps for each in- 
dividual CMC in Section p?. 1.1 1 and the Appendix. 

2.2. NRO CO( J =1-0) Data 

The CO (J =1 — 0) emission was observed in 
a 31' X 36' area (corresponding to 7.6 kpc x 8.8 
kpc) toward the disk of M33 using the NRO 45-m 
telescope ( Tosaki et al.|201l" ). The data reduction 
was also carried out with the NOSTAR reduction 
package. The angular resolution of the final map 
was 19'.'3 and the grid spacing 7'.'5. Further pro- 
cessing was done with MIRIAD ( |Sault etal.|1995| , 
in a manner similar to the CO(J = 3 — 2) map. 
We have convolved the CO(J =1 — 0) map to a 
common angular resolution of 25" and re-gridded 
it to 8'.'0 per pixel. The convolved map had an 
rms noise of 87 mK for the velocity resolution of 
2.5kms~^ We have also created CO(J =1 — 0) 
integrated intensity maps for each CMC (see Ap- 
pendix), integrating over the same velocity range 
as used for the individual C0( J = 3—2) integrated 
intensity maps. 



2.3. Optical and Infrared Data 

We have used the photometric optical dataset 
from the f/5V7 ground-based survey of local star- 
forming galaxies with the Kitt Peak National Ob- 
servatory 4 m telescope presented in |Massey et al 



( 2006 ) , and available on their ftp site ^ It contains 



a total of 146,622 stars in a field of 0.8 deg cen- 
tered on M33, which fully covers the regions that 
were observed in C0( J = 3-2) andCO(J = 1-0). 
The photometry in B VI bands is less than 1 % 
above 21.1 mag and even less than 10% at about 



23 mag. Refer to Massey et al. (20061 for further 



details about the observation, data reduction and 
photometric analysis. 

We have also compiled available Hii region cata- 



logs ( Courtes et al. 1987 Wyder et al. 1997 Hodge 



et al.||1999D. The Ha image from |Hoopes k, Wal 
terbos ( 2000 ) is used to check the accuracy of the 



positions and extents of the M33 individual Hii 
regions in these catalogs. Refer to|Hoopes fc Wal- 



terbos ( 2000 ) for more details on the observations 



and reduction process. 

Finally, we have used the 24 /zm Spitzer data 



presented in Onodera et al. (2010) to correct the 



extinction in the Ha emission. The final combined 
mosaic image is approximately 1.1 x 1.2 deg^, and 
both the point-spread-function FWHM and grid 
sizes are set to 5'.'7. For more details on the reduc- 



tion process, refer to Onodera et al. (2010, and ref 



erences therein). Note that Onodera et al. (2010) 
have performed a pixel-by-pixel analysis and have 
not applied a local background subtraction to cre- 
ate the extinction-corrected data by combining 
the Ha image and the 24 data. The lack 
of background subtraction may cause a system- 
atic shift in the extinction-corrected Ha luminos- 
ity (e.g. Kennicutt et al.|[2009 ). In this paper, in- 
stead, we present a new approach in which the Ha 
and 24 /im luminosities are measured with circular 
apertures of several 100 pc radius and with local 
background subtraction applied (see Section|3.3[). 
We use the catalog of 24 ^m sources in |Verley et al. 
( 2007 ) , which is the only reference listing 24 fim 



sources to date, for the sake of comparing the po- 
sitions of Hii regions. 



^ ftp: / / ftp.loweII.edu / pub / massey /Igsurvey/ 
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3. Results 

3.1. The CO( J = 3 2) GMC Catalog 

The GMC catalog is obtained from our CO (J = 
3 — 2) data using the CPROPS package described 
in Rosolowsky & Leroy ( 2006 1 . Briefly, we search 



for compact emission in adjacent pairs of chan- 
nels with Tmb > 4: arms, the so-called cores, and 
then for emission with T^b > 2 (Tims connected 
to such cores. CPROPS attempts to account for 
the amount of flux below the 2 (Trms cutoff by 
linearly extrapolating the emission profile to the 
OKkms"^ intensity level. Finally, the emission 
identified with CPROPS must have a minimum 
velocity width of 2 channels (5 km s^ ^ ) with a min- 
imum peak of Tnib > 4 (Trms, SLud an equivalent size 
of at least the spatial resolution of 25" (~100pc). 
Therefore, GMCs with a narrower velocity width 
than 2.5 km s^^ (if any) will not be identified by 
this analysis. 

As a result, we have identified a total of 
71 GMCs, whose properties are summarized in 
Table [2j The columns indicate the cloud ID, 
the cloud position (R.A., Deck), the velocity, 
the peak intensity, the FWHM of the velocity 
width (AVfwhm), the major/minor axes and 
orientation of the GMCs, the CO(J = 3-2) 
luminosity (^^0(3-2))' ^irial mass (Mvir) and 
cross identifications with the CO(J = 1 — 0) 
and CO(J = 2—1) GMC catalogs of previous 
papers (Rosolowsky et al. 2007[ [Gratier et al 
2012). The cloud position and velocity are de- 
rived from an intensity-weighted mean over all 
the pixels of a GMC. The major/minor diameters 
(^maj, ^min) without beam deconvolution are also 
indicated. The virial mass, calculated under the 
assumption that each cloud is spherical and can 
be parametrized by a density profile p cx r~^, 
is given by Myj,. = 189 AV^y^^y^Rdeconv Mq 
(Solomon et al. 19871, where i?deconv is the de- 



convolved radius given by deconvolving the beam 
size (^beam) from the GMC size, i.e., i?deconv = 

v/[^maj-CaJ^/^[^Ln-^La„.]^/^- Notcthat 

the virial masses of 27 GMCs are not shown be- 
cause their minor diameters are smaller than the 
beam size. 

Among the 71 GMCs in our catalog, only one 
GMC (GMC-50) is a new identification. The 70 
GMCs have been previously identified in lower-J 



CO transition catalogs (Rosolowsky et al. 2007 
Gratier et al.|2012| . GMC-50 is located outside of 



the observed field in Gratier et al. ( |2012 ) and has a 
marginal detection in the C0{ J = 1—0) integrated 
intensity map in Rosolowsky et al. (20071, but 



was not identified in their catalog. A total of 10 
GMCs in our catalog are composed of more than 
two smaller GMCs in previous catalogs, due to 
their higher resolution data (~50pc; Rosolowsky 



et al. 2007 Gratier et al. 2012h. On the other 



hand, a few single GMCs in the Rosolowsky et al 



(2007)'s catalog appear to be composed of more 
than two GMCs in our catalog. This is because 
of our 8-15 times better sensitivity data, since de- 
composition of neighboring clouds is done based 
on peak differences exceeding 2 (Trms (Rosolowsky 
fc Leroy|2006 ). The ID numbers in these cases are 



marked with an asterisk in Table [21 Since 6 GMCs 
(GMC-18, GMC-28, GMC-47, GMC-51, GMC-58, 
and GMC-71) are at the edge of the observed fields 
and their extents are uncertain, they are not in- 
cluded in further analysis and discussion. 

Histograms of the physical properties of 65 
GMCs are presented in Figure [2] The AFfwhm, 
and Mvir are in the range of 



dcconv ; ^00{3- 



■2)' 



1.2-17.3kms-\ 12-157pc, (0.8-25.8) x 10^ K km s~ 
pc^, and (0.2-89.2) xlO^M©, respectively A 
lower limit for the virial mass is 1.2 x 10^ M^q, as- 
suming an effective radius of ^50 pc and AVfwhm = 
2.5 km s^^, which are close to the instrumen- 
tal limits. Average values are AVfwhm ~ 
8.3kms"\ -Rdcconv < 73 pc, Lco(3-2) ~ 5.1 x 
lO^Kkms^ipc^, and Mvir < 7.7 x 10^ Mq. Note 
that we have used the non-deconvolved radius for 
the GMCs whose minor axes were smaller than the 
beam size, and thus the averaged i?doconv and Mvir 
should be regarded as upper limits. These aver- 
aged values are within the typical range as traced 
by CO(J 



2010) 



1 — 0) in LGGs (Fukui & Kawamura 



Compared to the high-resolution CO(J = 3 — 2) 
observations towards clumps in LMC, the i?dcconv 
and Mvir of M33 GMCs are larger than the typical 
values of LMC clumps, i.e. i?dcconv = l.l-12.4pc 
and Mvir = 4.6 x 10^-2.2 x 10^ Mq ([Minamidaiii 



et al. 2008). This is likely only due to our more 



limited spatial resolution, because the AVfwhm 
of M33 GMCs is similar to that of the clumps at 
high resolution, AVfwhm of ^7kms^^. This sug- 
gests that the M33 GMCs obtained in this work 
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are composed of a single or a few small clumps at a 
similar velocity, but not a large amount of clumps. 

3.1.1. CO(J ^3-2) and CO(J = 1 - Oj distri- 
bution and -R3-2/1-0 

The C0( J = 3-2) and C0( J =1-0) inte- 
grated intensity maps and also the CO intensity 
ratio (^3-2/1-0) map for each individual GMC, 
with a common 25" resolution, are shown in pan- 
els (a)-(c) of Figures 13 ~ 16 and Figures 18 - 78 



in the Appendix. The integrated intensity maps 
of CO(J = 3 — 2) emission have been created 
by using the channel maps where the emission 
is above 2(Tch- The integrated intensity maps of 
CO(J = 1 — 0) emission have been binned over 
the same velocity ranges as that of C0( J = 3 — 2) 
emission. 

We find a general trend that the distribution of 
CO(J = 3 — 2) emission and peaks are similar to 
that of C0( J =1-0) (e.g., GMC-3 in Figure 19, 
GMC-5 in Figure 21). However, in some cases 
28%), the CO(J = 3-2) emission for some GMCs 
shows more compact distribution (e.g., GMC-1 in 
Figure [15) GMC-8 in Figure [l4|. In other GMCs 
(~ 11 %), the C0( J = 3—2) emission peak is offset 
from the CO(J = 1 — 0) emission peak by over 



half HPBW (e.g., GMC-1 in Figure[15) GMC-7 in 
Figure 23). A further description of the spatial 
comparison between CO emission and massive SF 



sites is presented later in Section 3.5 

We have calculated i?3-2/i-o by dividing the 
C0( J =1-0) map by the CO(J = 3-2) map, 
after masking the region where the C0( J = 1 — 0) 
emission is below 3 a in its integrated intensity 
map. Note that each i?3-2/i-o map is masked 
except inside the boundary of the GMC, to avoid 
confusion with neighboring GMCs. 

The i?3_2/i-o values across GMCs vary greatly. 
Some GMCs are found to have a gradient of 
i?3_2/i-0j while others have a relatively constant 
ratio. The averaged i?3-2/i-o is measured inside 
the boundary of GMCs and is provided in Col- 
umn (9) of Table [2) The CO(J = 3-2) emission 
for four GMCs (GMC-49, 56, 57 and 70) have been 
detected with significant signal-to-noise ratio but 
the C0( J =1 — 0) emission has not been detected 
at the same position. In case that the i?3_2/i_o for 
these GMCs is uncertain, we calculate the lower 
limit of -^3-2/1-0 the CO(J = 3 — 2) emis- 



sion peak (see captions in each figure) . Figure [s] 
shows the histogram of i?3_2/i-o for the 65 GMCs. 
The averaged i?3_2/i_o ranges from 0.18 to 0.89, 
with a mean of 0.43. High ratios, i?3_2/i-o > 0.6, 
are found at four GMCs: GMC-1, GMC-7, GMC- 
34 and GMC-42, which are located around the 
GHRs (NGC 604, NGC 595) and the vicinity of 
the galaxy center. 

3.2. Identification of YSGs and their Ages 

The color magnitude diagram (CMD) is a pow- 
erful tool to provide insight into the age of the 
stellar component in a region of interest. Next we 
describe the method to identify YSGs as well as 
how to estimate their ages using the CMD. The re- 
sults are then used to investigate the spatial corre- 
lation between the YSGs and the GMCs identified 
in Section l3Al 



3.2.1. Extraction of Young Stars 

We have created the V magnitude - {B — V) 
CMD (i.e., Mv-{,B - V)o diagram) for the disk 
of M33 using the photometry catalog described 



in Section 2.3 To obtain intrinsic absolute mag- 
nitudes and then accurate ages, the photometric 
measurements must be corrected for two sources 
of reddening: the foreground extinction from our 
Galaxy and the internal reddening due to the in- 
terstellar matter in the disk of M33. We adopt 
the foreground galactic extinctions in the B and 
V bands of 0.181 and 0.139 mag, respectively 
dSchlegel et aI1[T998| . 

The reddening E{B — V) has been measured to- 
wards some individual clusters in M33, and ranges 
from 0.06 to 0.3 mag, with a typical value of 0.10 
(iChandar et all 119991 iPark & Leel 120071). Due 



to the different density structures of the GMCs 
it is necessary to consider a variable internal red- 
dening correction. This is especially important if 
the stars are embedded or lie behind the GMC, 
as the visual extinction is then expected to be 

The exact cor- 



Pineda et al. 


2010) 


given line oi 


sight 



tain, but we have used the available E{B — V) 
values obtained for some clusters, and the typ- 
ical value of E{B - V) = 0.10 ([Chandar et al. 



1999 Park & Lee 2007) for the remaining stellar 
The ratio of total to selective extinction 
V) 



groups 
Ry = Av/E{B 



3.1 in our Galaxy (Sav- 
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age & Mathis 1979) is used to obtain the visual 
extinction Ay- 

We have used the theoretical isochrones in 
the Padova stellar population synthesis models 
Girardi et al. 2010[ and ref- 



(Marigo et al. 2008 



erences therein) to estimate the age of the stars. 



The adopted distance to M33 (840 kpc; Freedman 



et al. 2001 ) corresponds to a distance modulus 
(m - M)o = 24.6. Note that the metalhcity (Z) 
varies from 0.6 Zq in the central region to QAZq 



at 5 kpc galaxy radius (Bresolin 2011 ). 



Figure [ijja) shows a Mv-{B — V)o diagram 
for the stars in M33 with five stellar isochrones 



from the Padova models (Marigo et al. 2008 Gi- 



rardi et al. 2010 and references therein), assum- 
ing Z = 0.5 Zq, ages of 6, 10, 20, 30 and 100 Myr, 
and an average internal extinction correction of 



Ay ~ 0.22 ( 


Massey et al. 


2007| Bastian et al. 


2007 


). Massey et al. 


( 


2006 


2007) have esti- 



~ 40 % of stars have {B - V)o be- 
1.0, but My in the range of 14.6- 



mated that 
tween 0.3 — 

19.6 are likely foreground objects (iBastian et al 



2007). Taking into account photometric errors 
and foreground galactic extinction, the criterion of 
{B — V)o < 0.3 is chosen to avoid the foreground 
stars. The faint limit of the V^-band magnitude 
of My < 21.5 is set to select only young massive 
stars (< 100 Myr). The minimum (initial) stellar 
masses at My ~ 21.5 for ages of 6, 10, 20, 30, 
100 Myr correspond to 14, 13, 10, 8, and 5Mq, 
respectively. The orange boundary in Figure |4]j a) 
indicates our selection criteria, {B — V)o < 0.3 and 
My < 21.5, an area that should represent young 
stars. 

The selected young stars are plotted in Fig- 
ure |4]jb), where the point size in the plot is propor- 
tional to the stellar brightness in y-band. These 
young stars (< 100 Myr) are found to be dis- 
tributed in a structure with two main arms and 
several weak multiple arms, which is in good agree- 



ment with Ivanov & Kunchev (1985). 



3.2.2. 



Young Stellar Number Density Map, Iden- 
tification of YSGs and their Ages 



We describe the method to search for YSGs and 
how to estimate their ages. The term "YSGs" 
refers to the concentration of young stars selected 
using the criteria in Section [3.2.1] These will in- 
clude both (young) star clusters as well as OB 



associations. Generally, a star cluster is defined 
to be a gravitationally bound system of several 
stars or more, whose concentration is larger than 
that of the surrounding stellar background, with 
a typical size of a few parsec to tens of parsecs 
(e.g., Hwang fc Lee||2008 and references therein). 
The definition of an OB association is usually a 
single, loosely bound or unbound concentration of 
early-type luminous stars, typically extending sev- 



eral tens to over hundred parsecs (e.g. Gouliermis 
et al. 2000 and references therein). The typical 
masses of young (< 10 Myr) clusters and associa- 
tions are 10-10"^ Mq and lO'^-lO^ Mq, respectivel y 
(e.g., [Hunter et al.|2003(|Kaleida fc Scowen|2010" |. 
Here we do not aim to distinguish between them 
because our interest is to investigate the collective 
properties of YSGs associated with GMCs. 

In previous studies, clusters or associations 
have been selected to be as stellar groups just 
by visual inspection, and therefore the selection 
may have been subjective. Recently, a more objec- 
tive approach has been carried out by [Gouliermis 
et al. (2000 20101), who have performed a stellar 



cluster identification in NGC 6822 and LMC using 
the stellar surface density map, and have shown 
the hierarchical structure of blue stellar clusters. 

To select YSGs in an objective manner for M33 
as well, we define them as an excess in the num- 
ber of young stars per unit of area (n*) following 
Gouliermis et al.[ ([2000| [20101) 's approach. We set 



the unit of area to the same scale as the pixel size 
of C0( J = 3 - 2) map, a 8" x 8" box (correspond- 
ing to ^ 30pc X 30 pc), which is comparable to 
the typical GMC size. The pixel size of the stel- 
lar number density map is set to half of the pixel 
size of the CO(J = 3 — 2) map, 15 pc. Figure [5] 
shows the number density map of young stars per 
each 30 pc x 30 pc area in M33 disk. It becomes 
more clear than in Figure |4]^b) that most of the 
regions with n* > 2 stars per unit area are dis- 
tributed along the spiral arms of M33. The max- 
imum value of stellar number density in M33 is 
found in the center of NGC 604, n* = 28 stars 
per area. Over-density regions {n* > 15 stars per 
area) coincide with prominent Hii regions such as 
NGC 604, NGC 595, NGC 592, IC 135, IC 139, 
IC 140, and the galaxy center, as labeled in Fig- 
ure [S] 

Next we compare the number density map with 
other wavelength data. Figure [6| shows the M33 
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disk in CO(J = 1-0), 24 /im, Ha, and FUV. 
The background-subtracted FUV map is a repro- 



duction from Verley et al. (20101. Interestingly, 



the distributions of the molecular gas traced by 
CO(J = 1 — 0) in Figure |6]ja) and the young stars 
in Figure [s] exhibit slightly different patterns be- 
tween the northern and the southern sides of the 
disk: molecular gas is more abundant in the north- 
ern side of the disk, while the opposite is true for 
young stars. Note that this asymmetrical pattern 
of the young stellar distributions is not because 
the photometric depth is inhomogeneous over the 
disk, as the completeness limits are similar (below 
V ~ 21.5 mag) both in the northern and southern 
disk. 

The distributions of 24 fim and Ha emissions 
show multiple spiral arms similarly to the den- 
sity map of young stars, but their spatial distri- 
bution patterns are not always similar: the 24 /im 
and Ha emissions are enhanced at GHRs such as 
NGC 604 and NGC 595 compared to other disk 
regions, while the surface density map of young 
stars is more pronounced along the spiral arms 
(Figure |6jD, c). The distribution of the young star 
number density is remarkably similar to that of 
the FUV emission in Figure |6][d), which traces the 
spiral structure and reproduces the north-south 
asymmetry. This is consistent with FUV being 
a diagnostic of the younger (30-100 Myr) stellar 
population ( [Verley et aL]|2009[ |2010 1 . 

In this work, we adopt 5 young stars per 30 pc x 
30 pc area as the minimum surface stellar num- 
ber required to be classified as a YSG, and then 
we estimate its age from the CMD, assuming that 
stars are coeval. First, we search a region with a 
peak of n* > 5 stars per area around each GMC, 
then determine the extent of the YSG at a level 
of n* = 2 stars per unit area. A Gaussian is then 
fitted to the radial profile of the young stellar dis- 
tribution down to n* > 2, and the extent of YSGs 
is set as the FWHM of that Gaussian fit assum- 
ing a circular shape (^d). Second, we have created 
the Mv-{B — V)o diagrams for the stars contained 
within the radius of rd. Finally we estimate their 
age using the Padova models. 

As an example. Figure ^b.) shows the num- 
ber density map of young stars, overlaid on the 
CO(J = 3-2) map of the NGC 604 region. 
Three GMCs are identified, GMG-34, GMC-1 
and GMC-27, from north to south in this re- 



gion. We find three peaks of n* > 5 stars per area 
around these GMCs, at position (aj20007 ^J20oo) = 
(l'>34'°31!4, 30°47'48'.'0), (l''34'"33M, 30°47'03'.'5) 
and (li^34™33!0, 30°46'23'.'8) from north to south, 
which we refer to as YSG-71, YSG-73 and YSG- 
72, respectively. The radii of these YSGs are 
measured to be r^i = 9.5" (39pc), 21.2" (86pc) 
and 10.4" (42 pc), respectively. 

Next, we derive ages for each YSG by fitting 
theoretical isochrones to the CMDs of resolved 
stars. Figure [8] shows the CMD of the three iden- 
tified YSGs above, YSG-71, YSG-72 and YSG-73. 
Dot symbols represent the stars within the radius 
of and the lines represent the Padova theo- 
retical isochrone tracks, assuming a metallicity of 
0.5 Zq. Solid lines represent the upper limit of 
the estimated ages, regarding the magnitude and 
color of the brightest main-sequence (MS) star as 
MS turnoffs for a young star group. In this case, 
the more massive stars are assumed to have al- 
ready evolved from the MS lines at {B — V)o ~ 
and less massive stars are still in the MS. Dashed 
lines represent the lower limit of the estimated 
ages, considering that all stars in the YSGs are 
still in the MS at {B — V)o ^ 0, without any more 
massive stars in the group evolving to red giant 
branches. Although chances are low that these 
massive stars are detected at redder colors, be- 
cause the more massive the stars, the faster they 
evolve, no detection of stars at {B — V)o > 
can suggest that more massive stars have not ex- 
isted before. Note that in either case, we have 
used only stars at {B — V)o < 0.3 to avoid fore- 



ground stars (Section 3.2.1). The measurable age 
limit of these YSGs using the Padova models is 
3 Myr. In this way, the estimated ages of the three 
YSGs YSG-71, YSG-72 and YSG-73 are in the 
range of 3-28 Myr, 3-28 Myr and 3-5 Myr, respec- 
tively. The first two YSGs include a smaller num- 
ber of stars (11 and 21 O stars for YSG-71 and 
YSG-72, respectively) and mostly less bright stars 
(My > 19 mag), which makes a larger uncertainty 
in the estimations. 

Figure [7][b) shows the distribution of the young 
stars, overlaid on the Ha map (grayscale). The 
size of the blue circle representing each star is pro- 
portional to the brightness in y-band, that is, the 
larger, the brighter. Many bright young stars are 
found in YSG-73 and concentrated in the center 
of the Hii region. In fact, NGC 604 is known to be 
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the most luminous and massive Hll region in M33, ues of OB associations (e.g., Bica &: Schmitt|1995 



containing more than 200 0-type stars ( 


Drissen 


et al.||1993||Hunter et al.||1996| |Gonzalez ] 


Delgado 


& Perez||2000 Bruhweiler et al.||2003). The other 



Battinelh et al. 19961 Kaleida & Scowen 20101 



two YSGs (YSG-71 and YSG-72) do not have cross 
identifications in previous cluster catalogs, but are 
included in previous star complex and OB associa- 
tion catalogs (No. 46 in Ivanov|2005 Humphreys 
[fc Sandagc 



3.2.3. The Catalog of YSGs 

We have identified 75 YSGs in M33 and have 
estimated their ages as described in Section [3. 2. 2[ 
Note that only YSGs associated with the identified 



GMCs are included (see Section 3.3). 

Figures [9}|ll] show the CMDs of 75 YSGs with 
the adopted Padova theoretical isochrone tracks. 
Similarly to Figure |8j solid and dashed lines rep- 
resent the upper and lower limits of the estimated 
ages. Blue dots represent the stars within the rd 
of the YSGs. In Figures [9l[l0l and[lTJ we show the 
CMDs of 18 YSGs younger than at most lOMyr, 
17 YSGs older than at least lOMyr, and the re- 
maining 40 YSGs, respectively. 

Table |3] shows the properties of the YSGs iden- 
tified in our catalog, including positions, size (rd), 
number of O stars, T^-magnitude of the brightest 
O star in the YSG, applied reddening correction, 
estimated age and cross identifications. The num- 
ber of O stars are counted only for those within 
our criterion on the My-(i? — V)q diagrams, i.e., 
{B - V)q < 0.3 and My < 21.5. Note that this 
stellar age provides just an upper limit for the 
YSGs, partly because reddening correction is un- 
certain. Twenty-three YSGs do not coincide with 
previously identified clusters, but are included in 



star complex and OB association catalogs ( Ivanov 
2005[ [Humphreys fc Sandage||1980[ ). The remain 
ing are found to have cross identifications with pre- 
vious cluster catalogs (SM; Park fc Lee]|2007 San 
Roman et a"L]|2010[ ). 

The radii of the YSGs range from 6'.'6 to 22'.'4 
(27-91 pc) and the average is 11" (46 pc). The 
number of O stars in a YSG ranges from 9 to 169 
stars. The total stellar mass of the YSGs are esti- 
mated to be lO^-s - IO^ '^Mq, assuming that the 
YSGs are characterized by a Salpeter IMF (slope 
-2.35), with masses spanning 0.1-58 M©. These 
sizes and masses are comparable to the typical val- 



It is possible to estimate the stellar spectral type 
of the primary ionizing star contained in a YSG 
from the brightest V-magnitude and {B — V)o 
color. The most massive stars are found in YSG- 
21 and YSG- 73 and their masses are estimated to 
be 51 Mq, corresponding to 03.5V (Martins et al 



2005D . The brightest stars in YSG-39 and YSG- 
60 are less massive than in any other YSGs, 8 Mq 
(corresponding to B2V). The latter is compara- 
ble to small clusters in our Galaxy and LMC that 
can be ionized by a single mid-0 or BO star (e.g., 
Panagia & Ranieri||1973| |Wilcots||1994|. The ages 



of the YSGs are estimated to be in the range be- 
tween 4 ± 1 Myr and 31 ± 19 Myr, with an average 
of 12 ± 5 Myr. 

We can compare the ages with other catalogs. 
Although the YSGs in our catalog do not fully 
overlap with the clusters whose ages have been 
previously determined (SM; PPL; SSGH), only 5 
YSGs coincide, YSG-1, YSG-13, YSG-41, YSG- 
63, and YSG-67. The previously derived ages are 
25t?!Myr (PPL), 13Myr (SM), 16^^ Myr (PPL) 



PPL) 



-14 ■ 

or 21 Myr (SM), 5 Myr (SM) or 13l5 Myr 
and 40I24 (PPL), respectively, while they are 
9 ± 2 Myr, 11 ± 5 Myr, 9 ± 2 Myr, 5 ± 2 Myr 
and 14 ± 11 Myr in our estimation. Compared 
to the results for YSG-1, YSG-41, YSG-63 and 
YSG-67 in PPL who have used a similar method, 
we notice that our estimations of ages are slightly 
younger than them but consistent within the un- 
certainties. This slight diff'erence is mostly be- 
cause they have used a small member star selection 
radius (within 1-2") of a compact stellar cluster 
and also because they have applied different the- 
oretical isochrones in the Padova models, corre- 
sponding to Z = 0.2 Zq. Besides, even though 
a cluster in their catalogs is associated with one 
of our YSGs, in general it only represents a small 
portion of the YSG, and thus it is not clear in 
principle whether they represent exactly the same 
entity. For the ages of the other two YSGs (clus- 
ters) in known GHRs (NGC 595 and NGC 604), 
YSG-2 and YSG-73, they had been previously es- 



timated to be 4- 6 Myr ( [Drissen et al.[1990[ [Malu 
muth et al.|1996 l and 3-5 Myr, respectively, using 
optical spectroscopy combined with an instanta- 
neous starburst model (e.g., Gonzalez Delgado &: 
Perez[[2000| . The ages of these young clusters are 
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in good agreement with our results, 6 ± 3 Myr and 
4 ± 1 Myr, respectively. 

3.3. Hii Regions and YSGs Associated 
with GMCs 

In this section, we investigate the associations 
of Ha, 24 /xm sources, and the identified YSGs, 
with the GMCs. These associations are listed in 
Table[4j Hii regions are tracers of the younger stel- 
lar population (< 10 Myr), while 24 /xm sources 
represent massive SF sites embedded in dust. 
When the extent of Hll regions, 24 /im sources 
and YSGs are confined within the boundary of a 
GMC, they are treated as being associated with 
the GMC. 



and the SFR (Calzetti et al. 2007). The L(Ha) 



The Hii region catalogs ( 


Boulesteix et al.|1974 


Courtes et al. 


1987 


Hodge et al. 


1999 


) include 



those Hii regions with Ha luminosities larger than 
1034.4 grgg-i ( |Wyder et aL]|l997| . Therefore, the 
sensitivity of the survey is high enough to detect 



the equivalent to a single 09V type star (Mar- 
tins et al.|20"05| , assuming a standard relation be- 
tween ionizing photon rates and Ha luminosities 



(Mayya & Prabhu 1996 and references therein). 



For the 24 /xm sources, the faint end of the 24 fj,m 



luminosity in the catalog (Verley et al. 2007) is 
^ lO"^^'^ ergs~^, corresponding just to an Hii re- 
gion illuminated by a single B1.5V star. Note 
that the 24 /xm source catalog includes discrete 
sources such as Hii regions, supernovae remnants 
and planetary nebulae, but most of the 24 /xm 



sources are classified as Hii regions (Verley et al. 



2007 ) , which contain embedded young stellar clus- 
ters with ages of 3-10 Myr ( Sharma et al.|[20TT ). 

The 24 /im sources have generally correspond- 
ing Ha emitting sources, but there are some cases 
of 24 /xm sources without corresponding Hii source 
(see the last column in Table [4]). All GMCs except 
GMC-60 are found to be associated with Hii re- 
gions and/or 24 /xm sources, while 51 GMCs are 
found to be associated with YSGs. In 10 GMCs 
no Hii region exists but 24 /im sources have been 
identified, while in 16 GMCs there are Hii regions 
without corresponding 24 /im sources. 

The Ha and 24 /im luminosities, _L(Ha) and 
1/(24 /im), and SFR for each GMC are also listed 
in Table 111 where the SFR for each GMC is calcu- 
lated from the L(Ha) and i(24 /im) using the rela- 
tion between the extinction-corrected Ha emission 



and L(24/xm) for each GMC, are measured with 
a circular aperture of size equal to that of the as- 
sociated GMC (100 to 400 pc diameter). These 
aperture photometry measurements include local 
background subtraction. Determination of the lo- 
cal background is done by fitting the emission over 
the aperture with Gaussian functions. These aper- 
ture sizes are similar or larger than a typical Hii 
region (Ipc - 100 pc) and thus they may poten- 
tially contain several Hii regions. A more care- 
ful local background subtraction for these Hii re- 
gions would just result in even larger differences 
in luminosities between a small Hii region and a 
large one because the local background contami- 
nation tends to be systematically large at the faint 
The L(Ha) and L(24/xm) 

to 7.67 X lO^^ergs" 
^37 f„ /I A/i ^ in40„^„„-i 



end (Liu et al 



2011j). 
from 2.49 x lO^^ 



and 



range 

from 1.67 x 10^^ to 4.04 x 10''" erg s , respec- 
tively. The most luminous source both in Ha and 
24 /xm emission is found at GMC-1, correspond- 
ing to NGC 604. The SFR of the M33 GMCs 
is characterized by a large scatter, ranging from 
8.89 X 10-5 to 0.216 Moyr-i. 

We estimate the averaged Ha attenuation 
over a GMC (i.e., a few 100 pc scale) using 



i(24/xm)/L(Ha) ratio as a proxy (Prescott et al. 
2007[ [Kennicutt et al.||2007[ ). The Ha attenuation 
is given by = 2.5 log[l-H0.038 L(24/im) /L(Ha) 



( [Calzetti et al.|2005[ [PTescott et al.|2007[ ). We list 



Afia in Table Wl The is generally low (less 
than 1 mag) and in good agreement with previous 

[2011] ). 



studies (Grossi et al. 2010 Sharma et al. 



The averaged A^a is 0.4 mag {Ay ~ 0.5) and is 
consistent with the adopted extinction for indi- 
vidual YSGs in Section [3.2.1[ Relatively highly 
obscured sources {Ana > 2 mag) are found in 
three GMCs: GMC-6, 15 and 68. The ages for 
the YSGs associated with these GMCs may be 
underestimated. Note that highly obscured re- 
gions {A^ia > 3 mag) have been found in other 
nuclei of galaxies (Prescott et al. 20071, but we 



can discard that possibility in M33. 

3.4. GMC Types 

We classify the GMCs into four types according 
to the age of the associated YSGs and Hii regions: 

Type A: GMCs are not associated with Hii re- 
gions nor YSGs; 

Type B: GMCs are associated with Hii regions. 
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but not with any YSG; 

Type C: GMCs are associated with Hii regions 
and young (< lOMyr) stellar groups; 

Type D: GMCs are associated with Hll regions 
and relatively old (10-30 Myr) stellar groups. 

Note that averaged values between the upper 
and lower limits of the stellar ages are used for 
this classification. If a GMC is associated with 
several YSGs of different ages, the classification is 
done based on the youngest stellar group. 

TableHhsts the classification of each M33 GMC 
into one of these four types. Out of the 65 GMCs, 
1 (1%), 13 (20%), 29 (45%), and 22 (34%) are 
found to be Types A, B, C, and D, respectively. 
Table [5] summarizes these statistics. The Type C 
GMCs are the majority among all GMC types. 

The estimated age of some YSGs may be un- 
certain. We select 44 GMCs with accurate age 
estimations for their associated YSGs in order to 
check how robust the relative percentage is among 
the different types. In this regard, only YSGs 
whose age upper limit is smaller than 10 Myr (see 



larger than in LMC (26 %) . Note that the clusters 
with ages of 10-30 Myr in LMC are likely far away 



Figure 
Figure 



9^ and lower limit larger than 10 Myr (see 
10| are used. With this condition we ob- 
tain that 1 (2%), 13 (30%), 19 (43%), and 11 
(25 %) are classified as Types A, B, C, and D, re- 
spectively. The 44 selected GMCs are indicated 
in Table ID We confirm the trend in the relative 
percentages, being Type C GMCs the majority 
among all GMC types. 

A caveat in the Type B and C classifications 
is that the percentages may be partly influenced 
by the instrumental sensitivity in [Massey et al.| 
(2006)'s survey, which may limit the detection of 
some YSGs. Given the instrumental sensitivity 



of Massey et al. (2006 1 's survey, some of Type B 



sources possibly contain YSGs that have not been 
detected. However, the sensitivity restricts the 
contained YSGs to be less massive than 13 Mq 



(-B1V star) (Section 3.2.1), and Type C (YSGs 



with ages < 10 Myr) are typically more massive 
than that. 

The three GMC classifications in the LMC, 



Types I, II and ill (Kawamura et al. 2009) corre- 



spond to Types A, B and C in our classification, re- 
spectively (see also Table [5]). In the LMC, Type ii 
(Type B) represents the largest number among all 
GMC types, while it is Type C in M33. The per- 
centage of Type C GMCs in M33 (40-45 %) is also 



from their parent GMCs (Yamaguchi et al. 2001 



Kawamura et al. 2009 ) and thus the relative per- 



centage is largely unknown. 

In addition, these discrepancies can be partly 
explained by the different spatial resolution and 
observed line, as the classification of LMC GMCs 
have been done with 40 pc resolution in CO (J = 
1-0) ( |Fukui et al.|2008| ). First, the GMCs in M33 
are typically double the size of those in LMC. If 
a M33 GMC is composed of two different types of 
GMCs, for example Types B and C, whose sizes 
are comparable to those in LMC, then the GMC 
would have been classified as Type C. Therefore, 
the classifications will be affected. Second, re- 
garding the different tracers used, the M33 GMCs 
identified from CO(J = 3 — 2) data represent 
the denser and/or warmer molecular gas which is 
more directly linked to massive SF, while the LMC 
GMCs identified from CO(J = 1-0) data trace 
the bulk of the molecular gas. In fact, in the case 
of M33 GMCs identified by using C0( J =1-0) 
data, more than two-thirds have associated Hii re- 



gions (Engargiola et al. 2003), which is a similar 
proportion in LMC ( Kawamura et al!]|2009 ). 



3.5. Physical Properties of GMC Types 

Figure [12] shows the distributions of the line 
width (AVfwhm), deconvolved radius (i?deconv), 
CO(J = 3 — 2) luminosity (-^co(3-2))' averaged 
line ratio (i?3-2/i-o) and SFR of the 65 GMCs. 
Each row represents each GMC type. The shaded 
areas represent the histograms for the 44 selected 
GMCs with more accurate classifications (Sec- 
tion Table [6] summarizes the mean and stan- 
dard deviation of these variables for each GMC 
type. The i?3_2/i-o distributions of Type C and 
D show a peak over 0.4, while those of Type B 
are always smaller than the average value of the 
former types (i.e., i?3-2/i-o< 0.3). The SFR dis- 
tributions of Type C and D show a peak over 
2 X lO~^M0yr~^, while those of the other three 
types are smaller. We have used Kolmogorov- 
Smirnov (K-S) tests to check whether the distri- 
butions for the different types differ. We cannot 
discard that the AVfwhm, -Rdoconv and -^co(3-2) 
distributions for these three types arise from the 
same distribution. However, the test indicates 
that i?3_2/i-o and SFR distributions of Type C 
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and D are different to those of Type B (p-value 
< 0.01). The selected GMCs with more accurate 
classification, as shown in the shaded histograms 
(Figure 12), indicate that the average values of 



all quantities for Type C are the largest, followed 
by Type D, Type B, and Type A, in this order. 
This suggests that the Type C GMCs are rela- 
tively large and show the most active SF among 
all types. Note that all quantities of the single 
Type A object are smaller than for the other three 
types. However, additional Type A sources need 
to be identified to properly characterize the prop- 
erties of this class of objects. 

In the following subsections, we present a close- 
up view of the CO(J = 3-2), CO(J =1-0), 
i?3_2/i-o Ha maps of the largest CMC in 
each CMC type, in order to illustrate the different 
properties among the four GMC types. Figures [T3| 
- [16] show for four GMCs representative of each 
type: (a) the C0( J = 3 — 2) integrated intensity 
map, (b) CO(J =1 — 0) map, (c) the i?3_2/i-o 
map, and (d) the Ha images. GMC-60 is chosen 
as a representative of Type A, GMC-8 of Type B, 
GMC-1 of Type C, and GMC- 16 of Type D. The 
positions of Hii regions, the 24 /im sources, and 
the surface density map of young stars are also 
marked in these figures. 

The C0( J = 3 — 2) emission peaks are in gen- 
eral slightly offset from the C0( J = 1—0) emission 
peak as we mentioned in Section |3.1.1[ but they 
are closer to the location of massive SF sites such 
as YSGs and Hii regions. These differences of spa- 
tial distributions necessarily results in a relatively 
high R3-2/1-0 area in a GMC around such massive 
SF sites, and the gradient of i?3_2/i_o that gradu- 
ally increase in the direction toward such massive 
SF sites. 

3.5.1. GMCs without Massive SF Regions ( Type A ) 

Figure [13] shows the images of GMC-60 of 
Type A. No large variations in i?3-2/i-o 0-28 
on average over the GMC) or bright Ha emission 
spots in this GMC are found. Besides it is a rel- 
atively small GMC compared to the average size 
and mass of all identified GMCs (Table [6]). Since 
only one GMC is classified as Type A, it is unclear 
whether the Type A source is representative of a 
class of similar sources or just a statistical outlier. 
Additional Type A sources need to be identified to 
properly characterize the properties of this class of 



objects. 

3.5.2. GMCs with Hii Regions hut not YSGs 
(Type B) 

One-fifth of the GMCs in our catalog are clas- 
sified as Type B. We cannot assure significant dif- 



ferences in AVpwHM, size and L', 



CO(3-2) 



compared 



to those of Type C and Type D (see Figure 12 1. 
However, the i?3_2/i_o is smaller than the latter 
two types, with i?3_2/i_o ^ 0.3 on average, and 
never exceeding 0.5. In addition, the SFR is also 
lower than the latter two types. We find that all 
Type B GMCs share a common trend, in the sense 
that higher i?3_2/i-o inside GMCs are preferen- 
tially found close to Hii regions and 24 /im sources, 
although the peak of i?3_2/i-o is not always coin- 
cident with the central position of Hii regions. No 
large spatial variation in i?3_2/i-o o'^Gr the GMC 
is found. 

The associated Hii regions and 24 /im sources 
show a relatively random distribution over the 
GMC. In addition, these sources are relatively 
small both in extent and in luminosity. In fact, 
we find that they do not usually exceed L(Ha) ~ 
lO^^ergs"-^ or L(24/im) ~ 10^*-^ergs^\ which 



corresponds to a single 07.5V star (Martins et al 



2005 ) , assuming the conversion factor from ioniz- 



ing photon rates and extinction-corrected Ha lu- 
minosities in Mayya & Prabhu (19961. This sug- 
gests that the associated Hii regions are ionized 
primarily by a star as massive as a 07.5V star, 
or later type of OB stars accompanied by several 
other less massive stars. 

As an example of Type B GMCs, we present a 
close-up view of GMC-8 in Figure [14] This GMC 
is known as the most massive GMC as traced by 
CO(J = 1- 0) ( jEngargiola et al.|2003| [Rosolowsky 
et al.|200"7 ), but it is just identified in our CO(J = 
3-2) catalog as the eighth brightest GMC. The 
CO(J = 3 — 2) emission shows a more compact 
distribution than that of CO (J =1 — 0) and the 
emission peaks coincide. The i?3„2/i-o slightly 
increases at the center of the GMC or close to Hii 
regions and 24 /im sources, but no large variation 
is found all over the GMC. 



3. 5. 3. GMCs with Hii Regions and YSGs (Type C) 

About half of the GMCs in our catalog are clas- 
sified as Type C. The AVfwhm7 size and ^^0(3-2) 
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of this type are similar to those in other types, 
but -R3-2/1-0 and SFR are the largest among all 
(Figure 12 and Table [6]). There are gradients 
of i?3_2/i_o around massive SF sites. Peaks of 
i?3_2/i-o do not coincide in general with the po- 
sitions of these SF sites, as also found in Type B 
GMCs. The average i?3-2/i-o of this type, ^0.5, 
is usually larger than that of Type B, and even 
exceeds over 1.0 at its peak. In addition, unlike 
Type B, we find that the associated Hii regions 
usually exceed the equivalent I/(Ha) of a single 
07.5 star. 

As an example of Type C GMCs, we show the 
maps for GMC-1 in Figure [l5| The CO(J = 
3 — 2) emission shows a more compact distribu- 
tion than that of CO(J =1-0). GMC-1 over- 
laps with two YSGs, YSG-73 (corresponding to 
the GHR NGC 604) and YSG-72. The emission 
peak of CO(J — 3 — 2) is offset from that of 
CO(J = 1 — 0) and is located northward, closer 
to the largest YSG, YSG-73, and the center of 
the GHR. A high i?3-2/i-o (> 0.8) area is found 
on the southern side of YSG-73, where compact 
Hii regions are located ( [Churchwell fc Goss|[T999| 



Tosakiet al.[2007l ). High i?3_2/i_o (> 0.8) prefer- 
entially around the GHRs suggests that these ar- 
eas may be excited by young massive stars ( |Tosaki| 
et al.||2007 ). A gradient of i?3_2/i-o extends from 
its peak to the south. Panel (d) shows that the 
concentrations of many bright (massive) stars are 
located at the heart of the GHR. Note that it is 
not certain whether the YSG-72 is physically as- 
sociated with this CMC and NGC 604 because 
there is no spatial correlation between the YSG, 
Hii region, and i?3-2/i-o- 

3.5.4. GMCs with Hu Regions and 10-30 Myr 
YSGs (Type D) 

About one-third of the GMCs in our catalog 
are classified as Type D. Generally, the physical 
properties of this type is similar to Type C (see 
Figure 12 and Table [6]). We also find that the 
associated Hii regions are slightly less bright than 
those in Type C, but larger than those in Type B. 

Figure [T6| shows the maps of GMC-16, a repre- 
sentative of Type D GMCs. The C0( J = 3-2) 
emission shows a similar distribution to CO(J = 
1 — 0), and both CO peaks coincide. Again, the 
peak of i?3_2/i_o (~ 0.4 — 0.6) is preferentially en- 
hanced in the vicinity of the associated massive SF 



sites but not always coincident with them. Such 
relatively high i?3_2/i_o area is elongated from 
north to south and there is a gradient of -^3-2/1-0 
from the -R3_2/i-o peak in the middle of the CMC 
to the both west and east edges of the GMC. 

4. Discussion 

4.1. Evolution of the GMCs 

In this section, we discuss based on our results 
the possible evolution of a GMC by focusing on 
the dense gas formation occurring around the first 
generated stars, which presumably leads to SF. 
We also estimate the typical lifetime of a GMC in 
M33. 

Since there is no correlation between Hii regions 
or YSGs, we infer that Type A is at an evolution- 
ary stage before massive SF. In the case of LMC 
with 40 pc resolution and using the C0( J = 1 — 0) 



line (Fukui et al. 2008), Kawamura et al. (2009) 



have shown that about a quarter of all GMCs are 
Type A. Possible reasons that so few GMCs in 
M33 are classified as Type A are spatial resolution 
and sensitivity limit and the choice of the molecu- 
lar gas tracer, as explained in Section [3^ Type B 
GMCs, which are associated with Hii regions but 
not with YSGs, are interpreted as being at an evo- 
lutionary stage just after the formation of massive 
stars, such as a single 07.5V star or later OB type 
stars, that are still partly embedded to be seen 
in the optical as YSGs. We suggest that Type C 
GMCs being associated both with Hii regions and 
very young stellar groups (less than 10 Myr) are 
at an evolutionary stage with active massive star 
formation, as massive as earlier types than 07. 5V 
stars. Type D GMCs associated with Hii regions 
and relatively old (10-30 Myr) stellar groups are 
considered to be at an evolutionary stage where 
they have been continuously forming massive stars 
for over at least 10 Myr. 

We find that -R3_2/i-o is enhanced around the 
massive SF sites for Types B, C and D GMCs 



(Section 3.5.2 3.5.4 ). In order to roughly quantify 
the offset between them, we plot in Figure [T7| the 
number of pixels above a certain i?3_2/i_o (0.1, 
0.3, 0.5, 0.7 and 0.9) as a function of the distance 
to its closest YSG. The average sizes of the 75 
YSGs are shown in the same figure, for reference. 
The dashed lines represent the frequency distri- 
bution expected if the same number of YSGs are 
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distributed at random in the observed area. It 
shows that the distance to the nearest YSG be- 
comes smaher as R^_2/i-o increases, but it must 
be noted that the i?3_2/i-o peaks and the nearest 
YSGs do not coincide in general. 

We can assume that i?3-2/i-o is a good tracer 
of density when i?3_2/i-o < 0.7, if kinetic tem- 
perature is low (^20 K), as explained in IMuraoka| 
et al. (2007). This suggests that dense molec- 



ular gas fraction is enhanced in the vicinity of 
previously formed YSGs. For regions with high 
i?3_2/i-o > 0-8, the molecular gas can have a 
higher kinetic temperature due to the nearest 
GHRs and YSGs (within ^50 pc). This is in good 
agreement with the suggestion by [Wilson et al.| 
( 1997 ) that the GHRs may raise the kinetic tem- 



perature of the neighboring (< 100 pc) molecular 
gas. Therefore, the enhancement of dense molecu- 
lar gas fraction around the previously formed mas- 
sive stars imply that once SF starts, dense gas for- 
mation is favored around the SF sites, which will 
then lead to the next generation of SF. 

Next we estimate a lower limit of the lifetime 
of a GMC from the age of the oldest stellar group 
contained within the GMC, assuming that the 
GMCs and YSGs are being formed in a nearly 
steady evolution. If the classification is done based 
on the oldest stellar group within a GMC, the 
frequency distribution of each type is 1 (2%), 13 
(20%), 19 (29%), and 32 (49%) out of 65 selected 
GMCs. The YSGs associated with Type C and 
Type D GMCs are 8 ± 4 Myr and 18 ± 10 Myr old 
on average, respectively. Therefore, we estimate 
that the timescale for each evolutionary stage is 
8 ± 4 and 10 ± 6 Myr in Type C and Type D, from 
the age difference between the associated YSGs. If 
we further assume that the timescale for each evo- 
lutionary stage is proportional to the number of 
GMCs, then Types A, B, C, and D are estimated 
to be ~1, 3-7, 5-10, and 8-17 Myr, respectively. 
As a result, the typical lifetime (from Type A to 
Type D) of a GMC with masses of > 10^ M© is 
roughly 20-40 Myr. Note that these results do not 
change significantly when we use only 44 GMCs 



with accurate ages for their YSGs (Section 3.4). 

In LMC, the lifetime of a GMC with mass as 
small as 5 X 10** Mq was estimated to be 20-30 Myr 



(Fukui et al 


|1999 


Yamaguchi et al. 


2001 


Kawa- 


mura et al. 


2009 


). The timescales of Types A, 



B and C in LMC (corresponding to Types i, ii 



and III) are 6 Myr, 13 Myr, and 7 Myr, respec- 
tively. The timescales for Type B in M33 are 
smaller than that in LMC, and for Type C they 
are similar or slightly larger. It is not appropriate 
to compare the lifetime for Type A in M33 and in 
LMC because of the limited number of GMCs in 
this type. No classification exists for LMC GMCs 
after Type C, thus we cannot compare (see Sec- 
tion 3.4 for a discussion). Taking into account 
that M33 GMCs are a few times larger than the 
LMC GMCs, M33 Type C and D GMCs are pos- 
sibly composed of smaller clouds. If that is the 
case, these small clouds in M33 Type C and D 
may not be coeval. This can explain the fact that 
the timescale of Type C GMCs in M33 are slightly 
larger than in LMC. 

As a conclusion, we propose a scenario as fol- 
lows. First, dense gas formation occurs at a cer- 
tain place in a GMC. After the molecular gas be- 
come dense enough to form stars (Type A), mas- 
sive SF occurs in such dense regions, forming Hii 
regions (Type B). About 4-8 Myr later, the first 
formed YSGs are not embedded anymore and be- 
come visible (Type C). Subsequent dense molecu- 
lar gas formation occurs due to the effect of previ- 
ously formed YSGs and the next generation stars 
are born in such a dense gas (Type D). In this way, 
the SF propagates from initially generated stars, 
until the reservoir of molecular gas is exhausted. 
Such a continuous SF can last at least 10-30 Myr 
in a GMC with a typical mass of > 10^ Mq dur- 
ing their lifetime of 20-40 Myr (after the dense gas 
formation) . 

5. Summary 

We present a GMC catalog using wide field 
(121arcmin^ in total) and high sensitivity {la — 
16-32 mK in Tn^b for a velocity resolution of 
2.5kms~^) CO(J = 3 — 2) maps of the nearby 
spiral galaxy M33 taken with the ASTE 10-m 
dish, and a complementary new catalog of YSGs 
for which we have estimated the ages. 

We summarize our main results as follows: 

1. We identify 71 CO(J = 3-2) GMCs from 
the C0( J = 3—2) data. We discard from the 
analysis 6 GMCs which are at the edges of 
the observed field. The remaining 65 GMCs 
are characterized by the deconvolved sizes in 
the range of 12-157 pc in radius and virial 
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masses of 1.5 x 10^-8.9 x 10^ Mq, which arc 
comparable to those of GMCs in our Galaxy. 

2. We identify 75 YSGs from the excess of the 
surface density map of yoimg stars that are 
associated with the identified GMCs. The 
total number of O stars in a YSG spans from 
9-169, corresponding to total stellar masses 
of lO''"'' 10*''^ Mq, assuming that they have 
a Salpeter IMF. The YSG ages are also es- 
timated using the Padova model, and are 
found to be in the range of 4- 31 Myr (12 Myr 
on average). These values are comparable to 
those of typical OB associations. 

3. We compare the GMCs with the distribu- 
tion of YSGs and classical Hii regions from 
a compilation in the literature. The identi- 
fied 65 GMCs are successfully classified into 
four categories: Type A showing no sign of 
massive SF, Type B being associated only 
with relatively small Hll regions, Type C 
with both Hll regions and young (< 10 Myr) 
stellar groups, and Type D with both Hii re- 
gions and relatively old (10-30 Myr) stellar 
groups. Out of the 65 GMCs, 1 (1%), 13 
(20%), 29 (45%), and 22 (34%) are found 
to be Types A, B, C, and D, respectively. 

4. Prom a comparison of the distributions of 

the CO(J = 3 - 2)/C0(J =1-0) inten- 
sity ratio, -R3-2/1-O) with the YSGs and Hll 
regions, we find that the i?3_2/i-o peaks 
within a CMC axe slightly offset from the 
SF sites, but preferentially located around 
them. 

5. We interpret the four types as representing 
an evolutionary sequence of the GMCs. As- 
suming that the number of the GMC types 
is proportional to the timescale of each evo- 
lutionary stage, we roughly estimate that 
they are ^1, 3-7, 5-10, and 8-17 Myr, re- 
spectively for Types A, B, C, and D. This 
yields that the lifetime of a GMC with mass 
of > 10^ Mq is 20-40 Myr, similar or slightly 
larger than the lifetime of LMC clouds. 

6. Finally, we propose a scenario where after 
molecular gas becomes dense enough to form 
stars, massive SF occurs in such dense part. 
Then the molecular gas around the first gen- 
eration stars becomes dense and forms the 



new generation of stars. In this way, the 
SF propagates from the initially generated 
massive stars and continues during the GMC 
lifetime. 
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RA (12000) 



Fig. 1. — The CO(J = 3 — 2) integrated intensity map with 25" resolution lOOpc), overlaid on the Ha 
image. Color range spans from 1 to llKkms^^ and contour levels are 1, 3, 5, 7, and 9Kkms^^. The 
rectangular boxes indicate the eight observed regions in CO(J — 3 — 2), and the number identifies each 
region in Table [l] The cross symbol in Box 1 represents the galaxy center. 
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Fig. 2. — Histograms of physical properties of GMCs as traced by C0( J = 3 — 2): (a) deconvolved radius 
(-Rdcconv), (b) Full Width Half Maximum (FWHM) hue width (AVfwhm), (c) luminosity (-^co(3-2))' ^"^^ 
(d) virial mass (Mvir). The shaded areas indicate the detection limit. The dashed lines represent the average 
values: i?doconv ^ 73 pc, AVpwHM — 8.3kms~\ ico(3-2) — 5.1 x 10^ Kkms~^ pc^, and Mvir ^ 7.6 x 10^ 

Mq. 
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Fig. 3. — Distribution of average i?3_2/i-o over each GMC. The GMCs at the edge of the observed fields or 
with underivable i?3_2/i_o were discarded. The dashed line represents the average value, 0.43. 
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Fig. 4. — (a) The Color-Magnitude Diagram of the entire disk of M33. The five lines represent isochrone 
tracks of 6Myr (blue), lOMyr (cyan), 20Myr (green), 30Myr (yellow) and 100 Myr (pink). A metallicity 
of 0.5 Zq was used. The orange boundary (horizontal and vertical lines) indicates our color and magnitude 
selection criteria, {B — V)o < 0.3 and My < 21.5. (b) The spatial distribution of identified young stars, 
selected using the same color and magnitude criteria. Note that the point sizes in this plot are proportional 
to the stellar brightness in V-band. 
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Fig. 5. — The number density (n*) map of yomig stars shown in Figure Qb), measured in 30 pc x 30 pc area. 
The color range spans from 2 to 26 stars per area, in 2 stars per area color steps. This corresponds to 0.002 
- 0.031 stars per pc^, in 0.002 stars per pc^ steps. The number density of rt* < 2 stars per area is not shown. 
Prominent Hii regions are also labeled for reference. A cross symbol in the middle of the image represents 
the galaxy center. 
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Fig. 6. — (a) C0(1 — 0) integrated intensity map (see Section 2.2). Contour levels are 1, 2, 3,..., 8 a, where 



the noise level la — l.TKkms ^. Images of (b) 24 /im, (c) Ha, and (d) FUV (See Sections 2.3 and 3.2.2 ) 



21 




■ 1 ■ ■ ■ ■ • : 

Ih34m39s 36s 33s 30s 27slh34m39s 36s 33s 30s 27s 
RA (J2000) RA (J2000) 



Fig. 7. — (a) The stellar density map toward NGC 604, overlaid with CO(J = 3 — 2) contours. Contour 
levels are 1, 2, 3, ... and 13Kkms"^. YSG-XX and GMC-XX are the IDs from our young stellar group 
and CO(J = 3 — 2) GMC catalogs. The crosses are the Hli regions from Courtes et al. ( |1987 ) (red), and 
Hodge et al. (1999) and Wyder et al. (19971 (orange). The green triangles are the 24 ^m sources in [Verley 



et al.| (2007). The color range is the same as in Figure [sj (b) The distribution of individual young stars, 



overlaid on the gray scale image of Ha emission of the NGC 604 region. The blue circles represent young 
stars which are selected with our color and magnitude criteria (see Figure [4| and their sizes are proportional 
to the stellar brightness in y-band. Note that only the stars inside the radius of 90" from the center of the 
NGC 604 are shown. 
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Fig. 8.— The color-magnitude diagrams of the three identified YSGs around NGC 604: YSG-71, YSG-72 
and YSG-73 (NGC 604), from left to right panel. The dots denote the stars within the selection radius, 9'.'5 
(39 pc), 10'.'4 (42 pc) and 21'.'2 (86 pc), respectively, which correspond to the extent at a n* = 2 stars per area 
level. The solid and dashed lines are the Padova theoretical isochrone tracks of the upper and lower limit of 
the estimated age, respectively. The identification number in the catalog and the derived age of each YSG 
are shown at the upper left and lower right corners of each panel, respectively. 
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Fig. 9. — The color-magnitude diagrams of 18 YSGs (among the 75 in the catalog), whose upper limit ages 
are estimated to be less than lOMyr. Dots denote the stars within the selection extent (rd) in Table [3j The 
lines and labels are the same as in Figure [8) 
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Fig. 10. — Same as Figure |9] but for 17 YSGs (among 75 in the catalog), whose ages are estimated to be 
older than lOMyr. 
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Fig. 11. — Same as Figure|9]but for the remaining 40 YSGs not included in Figures [9] and 10 A total of 25 
of these are estimated to be younger than 15 Myr. 
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Fig. 12. — Histograms of AVfwhm, ^dcconv, ^co(3-2)' -^3-2/1-0 SFR for the 65 GMCs. (a)-(d), (f)- 
(j), (k)-(o), and (p)-(t) panels show the properties of the GMC Type A, Type B, Type C, and Type D, 
respectively. The shaded areas indicate the detection hmits: AVfwhm — 2.5kms~^, i?deconv = 50 pc, and 
ico(3-2) = 10"* Kkms"^ pc^. Note that i?deconv is not derived for 31 GMCs because the size cannot be 
derived for those with minor axis less than the beam size, and thus we set the non-deconvolved radius as 
the upper limit. The (blue) filled histogram indicates the selected 44 GMCs which are associated with the 
YSGs with confirmed ages (see details in Section 3.4). 



26 




Fig. 13.— The distribution of (a) CO(J = 3-2) and (b) C0( J =1-0) emission for the single Type A 
GMC identified in M33, GMC-60, overlaid on the surface density map of young stars (same as Figure l5|. 
The integrated velocity ranges from Vlsr = —171 to — 161kms~^. Contour levels are 1, 2, 3, 20Kkms . 
Dashed contours indicate 0.5 Kkms~^ level. Note that the contour levels below 2Kkms~^ are not shown 
for the CO(J = 1 — 0) map. The red and orange cross signs, and green triangle signs are the same as in 
Figure [7j (c) The i?3_2/i_o map for GMC-60. Contour levels are 0.4, 0.6, 0.8, and 1.0. The color range 
spans from to 1.2, and is shown on the right side of the panel, (d) The distribution of Ha emission around 
GMC-60. The blue circles represent young stars which are selected using our criteria (Figure |4^). Their 
sizes are proportional to the stellar brightness in V-band. 



27 



o 
o 
o 



HI 

□ 



so'oo" ''(a) CO(3-2)J ^ v 

40" 



20" 



GMC-8 : Type B 



49'00" 



40" 



+30''48'20" 



50'00" 



o 
o 
o 

tN 



□ 



40" 



20" 



49'00" 



40" 



+30*'48'20" 




14s 12s 



10s 08s Ih34m06s 

RA (J2000) 



12s 10s 08s Ih34m06s 

RA (J2000) 



Fig. 14. — Same as Figure 13 but for GMC-8 representing Type B. The integrated velocity range used in 
the moment maps is T^lsr = —264 to — 239kms~"'^. 
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GMC-1 : Type C 
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Fig. 15. — Same as Figure 13, but for GMC-1 representing Type C. The integrated velocity range used in 
the moment maps is T^lsr = —259 to — 229kms~^. 
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GMC-16 : Type D 



SO'OO 



o 
o 
o 

u 

D 



49'00" 



+30M8'20" 



o 
o 
o 

(V 

o 
u 
01 

□ 



+30M8'20 




04s 02s 34m00s 58s Ih33m56s 

RA (J2000) 



02s 34m00s 58s Ih33m56s 

RA (J2000) 



Fig. 16. — Same as Figure 13 but for GMC-16 representing Type D. The integrated velocity range used in 
the moment maps is T^lsr = —264 to — 244kms~"'^. 
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450 




distance {pc) 

Fig. 17. — Frequency distribution of the projected distances from the nearest YSG of the pixels where 
i?3_2/i-o is above 0.1 (blue hue), 0.3 (red), 0.5 (purple) , 0.7 (black) and 0.9 (green). For comparison, the 
dashed lines shows the frequency distribution of the distance when the YSGs are distributed randomly. The 
vertical (black) dashed line indicates the averaged size of the YSGs, for reference. 
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Table 1 

Properties of the ASTE CO (J = 3-2) observations 



Region No. 


Map Center 


Area 


Int. Time 






R.A., 


Decl. 


[arcsec, arcsec] 


[hr] 


[mK] 




(J2000.0) 




(1) 


(2) 


(3) 


(4) 


(5) 


1 


l''33'°44-!4, 


30°38'12''2 


310 X 335 


22.5 


16 


2 


l''34'°08"5, 


30°37'36''4 


263 X 240 


8.7 


19 


3 


l*" 34" 12=8, 


30° 34' 28'.' 9 


150 X 135 


1.3 


25 


4 


l*" 33" 38! 5, 


30°32'45'.'0 


200 X 440 


21.8 


18 


5 


l''33"35!0, 


30°41'48'.'0 


120 X 120 


1.9 


32 


6 


l''34"05!5, 


30°48'02''9 


250 X 250 


12.8 


19 


7 


l''34"33!2, 


30°47'06'.'0 


210 X 210 


10.5 


25 


8 


l''34"38?5, 


30°40'55'.'0 


180 X 150 


7.0 


18 



Note. — Col. (1): No. of the eight observed regions in order of increasing galaxy radius. 
Col. (2): Central position of each region. Col. (3): Size of the observed region. Col. (4): 
Integration time for each region. Col. (5): r.m.s. for a 2.5 kms~^ channel, measured in an 
emission-free region. 
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Table 3 

m33 young stellar group catalog 



YSG ID 


R.A., 


DecL 




rci 


Nostar 




E(B - V) 


Age 


Cross ID 




(J2000.0) 




( ) 




(mag) 


(Myr) 






(2) 






(.4) 


\°) 


(,6J 


(I) 


\P) 


1 


0l''33""30!2, 


30°31'44 


'9 


8.3 


21 


17.62 


CIO"" 


7-11 


SM464, PL12, SSA1272 


2 


0l''33'"33!6, 


30°41'27 


'6 


12.5 


67 


16.01 


0.20" 


3-8 


NGC595, PL57 


3 


0l''33"'34=;3, 


30°32'05 


'8 


12.6 


30 


17.94 


0.10 


6-12 


SSA1261, SSA1267, SSA1272 


4 


0l''33'"35!0, 


30°37'02 


'4 


7.4 


17 


18.84 


0.10 


17 - 28 


SSA1284, SM113 


5 


0l''33'"35';2, 


30°39'18 


'1 


12.5 


18 


18.54 


0.10 


11 - 17 




6 


0l''33'"35!7, 


30°42'32 


'2 


8.9 


18 


18.96 


0.10 


14 - 25 


SSA1285 


7 


0l''33'"35'!9, 


30°36'29 


'0 


9.8 


24 


18.38 


0.10 


3-17 


SM128 


8 


0l''33"'36!4, 


30°31'49 


'1 


8.0 


15 


19.70 


0.10 


22 - 39 


SSA1294 


9 


0l''33'"36!8, 


30°36'35 


'2 


8.0 


14 


18.05 


0.10 


3-15 




10 


0l''33"'38!0, 


30°32'01 


'7 


< 10.0 


18 


19.16 


0.10 


3 - 35 




11 


0l''33'"38"7, 


30°32'07 


'9 


8.0 


11 


18.60 


0.10 


6-19 




12 


0l''33'"39"3, 


30°38'04 


'0 


9.5 


42 


18.28 


0.10 


3-17 


SSA1341, SSA1350 


13 


0l''33'"39!6, 


30°32'37 


'2 


17.2 


67 


18.11 


0.10" 


6-15 


SSA1343, SSA1345, PL139 


14 


0l''33"'40!5, 


30°38'28 


'0 


8.7 


28 


16.52 


0.10 


5-8 


SSA1344 (SM545), SSA1361 




















(SM154) 


15 


0l''33'"42!2, 


30°40'21 


'5 


12.1 


29 


19.36 


0.10 


22 - 28 


SSA1394, SSA1407 


16 


0l''33'"42!3, 


30°33'00 


'2 


9.3 


22 


17.64 


0.10 


7-11 


SSA1406 


17 


0l''33'"43';i, 


30°32'47 


'6 


11.9 


42 


17.64 


0.10 


7-11 


SSA1415 


18 


0l''33"43!5, 


30°39'05 


'7 


17.0 


65 


16.40 


0.10 


3-7 


SSA1430, SSA1416 


19 


0l''33"'44!3, 


30°31'48 


'0 


13.8 


44 


16.87 


0.10 


3-8 




20 


0l''33"'44!8, 


30°33'08 


'2 


11.4 


23 


18.94 


0.10 


19 - 35 




21 


0l''33'"44!9, 


30°36'31 


'l 


16.4 


63 


16.46 


0.10 


3-7 


SM565 


22 


0l''33"'44!9, 


30°35'53 


'5 


9.9 


25 


16.78 


0.10 


4-8 




23 


0l''33'"46!3, 


30°36'52 


'0 


14.4 


42 


15.64 


0.10 


5-6 


SSA1472 (SM567) 


24 


0l''33"46!7, 


30°36'26 


'9 


14.4 


47 


18.09 


0.10 


10 - 17 




25 


0l''33"M7!5, 


30°39'41 


'8 


11.3 


31 


18.99 


0.10 


14 - 22 


SSA1493 


26 


0l''33"M7^5, 


30°38'36 


'4 


9.0 


30 


18.23 


0.10 


7-15 


SSA1495 


27 


0l''33'"48!l, 


30°33'02 


'3 


15.1 


51 


17.42 


0.10 


6-12 


SSA1499 


28 


0l''33'"48!7, 


30°38'45 


'0 


8.0 


17 


19.25 


0.10 


11 - 35 


SSA1515 


29 


0l''33"'48!7, 


30°39'30 


'0 


14.0 


53 


18.06 


0.10 


5-17 


SSA1505 (SM572) 


30 


0l''33'"49!l, 


30°39'45 


'3 


13.4 


59 


16.30 


0.10 


7-8 


SSA1516 


31 


0l''33'"49?3, 


30°38'22 


'5 


12.0 


49 


17.81 


0.10 


8-14 




32 


0l''33'"50!l, 


30° 37' 29 


'6 


9.4 


10 


18.42 


0.10 


3-22 


SSA1527 


33 


0l''33"'50!5, 


30°33'46 


'1 


10.0 


12 


18.95 


0.10 


15 - 25 


SSA1536 


34 


0l''33""51?2, 


30°39'55 


'8 


12.5 


73 


16.65 


0.10 


5-7 


SSA1550 


35 


0l''33"'51!6, 


30°37'22 


'5 


10.0 


20 


19.79 


0.10 


3-35 




36 


0l''33""51!6, 


30''38'51 


'0 


12.2 


84 


17.15 


0.10 


5-10 


SSA1554, SSA1546 


37 


0l''33'"52!3, 


30° 37' 56 


'0 


15.1 


39 


17.93 


0.10 


8-12 




38 


0l''33'"52?6, 


30°39'20 


'3 


10.0 


46 


16.42 


0.10 


3-7 




39 


0l''33"'52?8, 


30° 37' 15 


'0 


8.6 


9 


18.14 


0.10 


3-15 




40 


0l''33"'53!0, 


30°38'59 


'4 


9.5 


60 


16.91 


0.10 


3-8 




41 


0l''33""54!5, 


30°33'04 


'3 


12.2 


32 


17.17 


0.10" 


7-11 


SSA1586 (SM215, PL167), 




















SSA1598, SSA1585 


42 


0l''33""55!6, 


30°37'40 


'0 


10.0 


19 


18.81 


0.10 


12 - 25 


SM223 


43 


0l''33"'56!7, 


30°40'20 


'0 


10.5 


15 


19.36 


0.01 


6-35 


SM236 


44 


0l''33""57!4, 


30°35'30 


'0 


15.0 


41 


17.68 


0.10 


12 - 14 


SSA1636 (SM241) 


45 


0l''33"'57?9, 


30°48'46 


'2 


8.9 


11 


18.98 


0.15" 


7-25 


PL81 


46 


0l''33'"58!6, 


30°35'22 


'5 


11.0 


41 


16.41 


0.10 


3-7 


SSA1647 (SM476) 


47 


0l''33'°58!6, 


30°48'37 


'5 


10.0 


13 


18.74 


0.10 


5-22 


SM253 


48 


0l''33"59!0, 


30°49'09 


'2 


8.0 


10 


18.83 


0.15" 


14 - 28 


SM251 (PL80) 


49 


0l''33'°59!l, 


30°35'37 


'5 


10.0 


32 


17.51 


0.10 


5-12 




50 


0l''33'"59!l, 


30°35'47 


'2 


10.0 


28 


19.13 


0.10 


15 - 19 


SSA1654 


51 


0l''33'"59=2, 


30°40'07 


'5 


20.0 


99 


18.00 


0.10 


6-15 


SSA1652 


52 


0l''34"0l!5, 


30°36'20 


'6 


10.8 


38 


17.61 


0.10 


6-12 


SSA1693 


53 


0l''34"01'!6, 


30°37'17 


'0 


8.8 


20 


18.05 


0.10 


3-15 


SSA1700, SSA1682 


54 


0l''34'"0l!8, 


30°46'28 


'4 


16.4 


31 


18.10 


0.10 


10 - 17 


SSA1694, SSA1685 


55 


0l''34'"02!2, 


30°35'56 


'o 


13.2 


13 


19.21 


0.10 


3-15 




56 


0l''34'"02!3, 


30°38'40 


'6 


22.4 


80 


17.03 


0.10 


5-7 


SSA1711 (SM276) 


57 


0l''34""02?6, 


30°37'04('5 


8.3 


9 


18.38 


0.10 


3-15 
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Table 3 — Continued 



YSG ID 


R.A., Decl. 










E{B - V) 


Age 


Cross ID 






(J2000.0) 




(") 




(mag) 




(Myr) 




(1) 




(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


58 


0l''34' 


"03!2, 30°46'36' 


'7 


9.2 


9 


18.50 


0.10 


7-19 


SSA1715 (SM280) 


59 


0l''34' 


"03!3, 30°36'43'. 


'6 


9.5 


19 


18.13 


0.10 


3-15 




60 


0l''34' 


''05!7, 30°47'49'. 


'8 


10.7 


21 


17.77 


0.20" 


3-14 


SM297 (PL197) 


61 


0l''34' 


"06!9, 30°47'22'. 


'6 


13.2 


30 


16.44 


0.20" 


3-8 


SSA1748 (SM305, PL201) 


62 


01 ''34' 


''07!6, 30°39'24'. 


'4 


16.4 


69 


18.44 


0.10 


14 - 25 




63 


Ol''34' 


"09!3, 30°39'07'. 


'7 


14.9 


68 


16.16 


0.10" 


3-6 


SSA1791, SSA1795, SSA1777 
(SM321, PL98), PL97 


64 


0l''34' 


''ll!5, 30°36'12' 


'2 


17.4 


31 


17.75 


0.10 


8-14 


SSA1825 


65 


01 ''34' 


"14=0, 30°33'39'. 


'7 


7.2 


16 


17.95 


0.10 


5-15 


SSA1862 


66 


0l''34' 


"14!7, 30°33'58'. 


'5 


10.0 


32 


17.34 


0.10 


8-12 




67 


01 ''34' 


"15=1, 30°33'48^ 


'1 


7.5 


22 


18.90 


0.10" 


3-25 


SM478 (PL27) 


68 


0l''34' 


"15!2, 30°35'00' 


'0 


8.6 


14 


17.89 


0.10 


5-15 


SSA1894, SSA1895 


69 


01 ''34' 


"16=0, 30°33'37'. 


'5 


10.0 


29 


16.66 


0.10 


3-7 


SSA1906 (SM480) 


70 


0l''34' 


"16=8, 30°39'18'. 


'o 


8.3 


10 


19.76 


0.10 


12 - 50 




71 


0l''34' 


"3l!4, 30°47'48'. 


'o 


9.5 


11 


19.12 


0.10 


3-28 




72 


0l''34' 


"33=0, 30°46'23' 


'8 


10.4 


21 


19.07 


0.10 


3-28 




73 


0l''34' 


"33!4, 30°47'03'. 


'5 


21.2 


169 


15.90 


0.10 


3-5 


NGC604 


74 


0l''34' 


"39!5, 30°4l'37'. 


'1 


8.7 


16 


18.73 


0.10 


7-22 


SSA2167, SSA2168 


75 


0l''34' 


"40!1, 30°41'56' 


'0 


9.2 


19 


17.69 


0.10 


3-14 


SSA2169, SSA2164 



Note. — Col. (1): Young stellar group (YSG) ID, in order by right ascension . Col. (2): Central position R.A., Deel. of the YSG. 
Col. (3): Size of the stellar group, whieh approximately eorresponds to the boundary of n* — 2 stars per area (soc Scction [33^ . Col. 
(4): Number of O stars (more massive than 09V type star). Col. (5): V magnitude of the brightest member in the YSG. Col. (6): 
Applied red dening correction for each YS G. Col. (7): Estimated age of th e YSG. Col. ( 8): Cross identifications with other catalogs. 
SM-ra: from |Sarajedini fc Mancone| (2007) . PL-n: from [Park fc Lee|f2007[ . SSA-n: from |San Roman etlil] l |2010( . NGC-n: giant Hll 
regions. 

"From |Chandar et al!] l |1999[ , [Park &: Lee| ]2007[ , [Gonzalez Delgado &: Perez| | |2000[ l and jChurchwell fc Goss] l |1999[ . 
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Table 4 

Properties of M33 GMCs and massive star formation 



CMC 


Type 


YSG 


L(Ha) 


L(24 /im) 


A(Ha) 


SFR 




Hll regions / 24 /im source 


ID 




ID 


(erg s ) 


(erg s ) 


(mag) 


(Mg yr" 


) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


1 


C" 


YSG-72, 73 


7.67x10^® 


4.04 xio"" 


0.20 


2.16x10" 


-1 


NGC 604 


2 


D 


YSG-45, 47 


4.16x10^^ 


6.71x10^* 


0.52 


3.58x10" 


-3 


B-005 (V-426), B-074, C-Z233A 


3 


C" 


YSG-65, 66, 67, 69 


2.59x10^* 


1.96x10'*^ 


0.27 


1.05x10" 


-2 


C-Z258 (V-115), B-066 


4 


D 


YSG-64 


3.80x10'** 


1.63x10'*^ 


0.16 


8.73x10" 


-3 


B-709 (V-174), C-Z407 (V-195) 


5 




YSG-38, 40 


3.48x10^** 


1.51x10'*^ 


0.17 


8.09x10" 


-3 


C-Z160A (V-84) 


6 




YSG-5 


1.40x10^'^ 


2.13x10'*^ 


2.08 


1.13x10" 


-2 


B-033, B-628, B-203 (V-280, V-273), H-634 (V-234) 


7 




YSG-2 


2.89x10^° 


1.42x10*" 


0.19 


7.60x10" 


-2 


NGC 595 (V-184) 


8 






2.17X 10^'^ 


9.04x10'*^ 


1.03 


4.81x10" 


-4 


C-Z322. C-Z324, B-647 (V-434), B-302 (V507), V-407 


9 


C 


YSG-7, 9 


8.54x10^'' 


4.76x10'*** 


0.21 


2.55x10" 


-3 


C-Z112, C-Z214, V-168 


10 


C'' 


YSG-22 


6.94X 10^'^ 


4.07x10'*** 


0.22 


2.18x10 


-3 


C-Z219 (V-65). B-641, B-633 


11 


D'' 


YSG-33 


1.16X 10''** 


1.47x10'*** 


0.05 


8.01x10" 


-4 


B-047, B-014 (V-131) 


12 


C>' 


YSG-62, 63 


5.06X 10'" 


3.94x10^* 


0.28 


2.10x10 


-3 


B-6S9 (V-266), B-702 (V-240), C-Z251 (V-271) 


13 


C 


YSG-1 


2.47X 10^^ 


2.13x10^" 


0.31 


1.14x10 


-2 


B-211, C-Z109, C-Z045 (V-70) 


14 


C'' 


YSG-60, 61 


5.83X 10^* 


2.01 xio"^ 


0.13 


1.08x10" 


-2 


B-690, C-Z183, C-Z249 


15 


C*" 


YSG-18 


l.lOx 10^® 


1.89x10^* 


2.19 


1.01x10 


-3 


C-ZOlOA, B-044, B-03S (V-210) 


16 


D*" 


YSG-45, 48 


4.04X 10^'' 


1.66x10^* 


0.16 


8.90X 10" 


-4 


B-703 (V-431), B-669 (V-427), B-696 (V-429), C-224 


17 




YSG-56 


2.49X 10'^^ 


1.62xl0'*** 


0.24 


8.66x10 


-4 


C-Z405 (V-198) 


18 


c 


YSG-44. 46, 49, 50 


2.19X 10'^** 


2.55x10'*" 


0.40 


1.36x10 


-2 


C-ZOOIA (V-165), V-199 


19 




YSG-4 


1.95X 10'"^^ 


1.05x10'*" 


0.20 


5.62x10" 


-3 


B-251 (V-193) 


20 


C'' 


YSG-56 


3.80X 10'"^^ 


1.63x10'*" 


0.16 


8.73x10" 


-3 


C-Z405 (V-198) 


21 


D*" 


YSG-42 


1.47x 10'*^ 


1.89 x lo'*** 


0.43 


1.01 X 10 


-3 


IC142 


22 


D*" 


YSG-70 


1.61X 10'*^ 


3 83x1 o'*** 


0.70 


2.04X 10 


-3 


C-Z274 


23 




YSG 91 '?'^ 94 


2 19x 10'^^ 


1 14x 10^^ 


20 


6 1 Ox 1 n 


-3 


J_J — ^ ^ i_j — yfyjyj tV— XOU/, V— 


24 


Q 


YSG-27 


1.99X 10'*^ 


1.07x 10^^ 


0.20 


5 73 X 1 D" 


-3 


B-631 (V-127), B-019, V-113 


25 


c 


YSG-52. 53, 57, 59 


7.26X 10^^ 


6.94X 10^^ 


0.04 


3.81 X 10" 


-4 


B-672, C-Z245 


26 


Y) 


YSG-68 


9 78x 1 n^^ 


6.72 X 10^^ 


0.25 


3 59 X 1 D" 


-3 


C-Z322A {V-47). B-660, C-Z172 


27 


13 


YSG-72 


7.27x 10'^^ 


1.64x 10'^^ 


0.09 


8.84x 10" 


-3 


NGC 604 


28 




VSG-43 51 


S.lOx 10"^'^ 


9.60x 10'^"^ 


0.84 


5.11 X 10" 


-4 


H-744. H-731 H-740 


29 




YSG-30. 34 


1.23X 10'^^ 


3.42x10'^^ 


0.11 


1.84X 10" 


-3 


C-Z150, B-064, C-Z144 (V-72) 


30 


c 


YSG-3 


1.46X 10^^ 


1.19X 10^^ 


0.29 


6.35X 10" 


-3 


C-Z210 (V-60) 


31 






1 ^5 >c 1 0^^ 


7 (T^if 1 n^® 

1 • A\ X U 


0.12 


^ 7Sy 1 f)" 


-3 


R-m 7 ('V-23Q^ 


32 


D 




1.07x 10^* 


s Qi x- 1 n^® 

O • i7 X A. X U 


0.30 


4.76x 10" 


-3 


R-n2n {'v-'20'\^ R-fi37 


33 




YSG-18 


1.54x 10^* 


7 ^5 1 

1 ■ o o A, X u 


0.18 


^ Q4y 1 ft" 


-3 


±j \jtj\j I V ^ X u J 


34 


D 


YSG-71 


1.04x 10^® 


1 1 5 1 

• XO ^ X U 


0.08 


1.16x 10" 


-2 


NGC 604 B-652 


35 


B*" 




6.04X 10^^ 


4.35x 10^^ 


0.03 


2.41x10" 


-4 


C-Z347 C-Z353 


36 


D 


YSG-54, 58 


6.17x10^'' 


2.22x10^* 


0.14 


1.19x10" 


-3 


C-Z291, B-081 


37 


c 


YSG-16, 17, 20 


1.31x10^* 


5.76x 10^* 


0.17 


3.09X 10" 


-3 


B-061, C-Z139, C-ZOlO 


38 


D 


YSG-12 


2.04x10^® 


9.55x10^* 


0.18 


5.11x10" 


-3 


B-612, B-058 (■V-144), C-Z135, C-Z081 


39 


C 


YSG-3 


2.77x10^^ 


1.19x10^^ 


1.05 


6.34x10" 


-3 


C-Z210 ("^^-60), B-034 


40 


D 


YSG-8, 10, 11, 13 


2.07x10^'' 


2.56x10"* 


0.42 


1.37x10" 


-3 


B-200, B-206 


41 


C" 


YSG-36, 38, 40 


2.18x10^* 


9.18x10'*'' 


0.02 


5.24x10" 


-4 


C-Z160A (V-84), B-OlO 


42 


c" 


YSG-25, 29, 30 


9.66x10^'' 


1.29x10'*" 


0.45 


6.88x10" 


-3 


C-Z144 (V-72), B-037, C-Z221, B-066 


43 


c 


YSG-41 


3.11x10^® 


1.24x10'*" 


0.15 


6.65x10" 


-3 


B-092 (V-97) 


44 


c 


YSG-39 


6.61x10^^ 


3.35x10'*'' 


0.19 


1.79x10" 


-4 


C-Z227 


45 


D 


YSG-35, 37 


1.01x10^'' 


7.23x10'*'' 


0.26 


3.86x10" 


-4 


B-1501 


46 




YSG-15 


3.28x10^'' 


7.27x10'*'' 


0.66 


3.87x10" 


-4 


H-601, H-598 (V-283) 


47 




YSG-55 


2.19x10'** 


7.03x10'*'' 


0.01 


4.10x10" 


-4 


C-Z246A (V-122) 


48 


B" 




2.77x10^*" 


2.59x10"'' 


0.33 


1.38x10" 


-4 


H-1145 (V-277), H-1146, H-1266 (V-281) 


49 


D 


YSG-13 


1.21x10^** 


3.54x10'*** 


0.11 


1.90x10" 


-3 


B-050 (V-91) 


50 




YSG-6 


5.80x10^*^ 


1.36x10'*** 


0.09 


7.33x10" 


-4 


NGC 595 


51 


ga 




6.06X10'*'' 


7.85x10'*' 


0.05 


4.28x10" 


-4 


B-209, B-610, V-182 


52 


C" 


YSG-62, 63 


5.73x10'" 


1.54x10'*** 


0.11 


8.29x10 


-4 


C-Z182, B-689 (V-266) 


53 


B'' 




3.89x10^*^ 


2.63x10'*** 


0.25 


1.41x10" 


-3 


C-Z041 (V-69) 


54 


B" 




4.12x10'*'' 


1.48x10'*** 


0.14 


7.94x10" 


-4 


B-665 (V-438), V-398, V-407 


55 


B» 




9.58x10'*^ 


6.81x10'*' 


0.03 


3.78x10" 


-4 


B-302 (V-507), V-420 


56 


ci- 


YSG-18 


1.54x10"** 


7.35x10'*** 


0.18 


3.94x10 


-3 


B-038 (V-210), B-051 


57 


B^" 




8.55x10'*® 


3.89x10'*' 


0.17 


2.08x10" 


-4 


C-Z212 (V-422), H-769 (V-385) 
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Table 4 — Continued 



GMC 


Type 


YSG 


L(Hq) 


L(24 ^lm) 


A(Ha) 


SFR 




Hll regions / 24 /im sour 


ID 




ID 


(erg 




(crgs-i) 


(mag) 


(Mq yr- 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


58 


B" 




6.26 




2.50x10'" 


0.15 


1.34; 


<10 


-4 


C-Z146 (V-209) 


59 


D'=' 


YSG-20 


5.53 


<10^^ 


1.26x10^** 


0.09 


6.79; 


<10 


-4 


B-021, C-Z220, B-636 


60 


A" 




3.81 


<10^^ 


1.67x10^'' 


1.06 


8.89; 


<10 






61 


C" 


YSG-12, 14 


1.81 


<10^** 


9.55x10^** 


0.20 


5.11; 


<10 


-3 


B-612, B-058 (V-144) 


62 


B" 




1.09 


<10^^ 


8.03x10^'' 


0.27 


4.29; 


<10 


-4 


C-Z299 (V-402) 


63 


C 


YSG-74, 75 


2.68 


<10^** 


9.17x10^^ 


0.13 


4.92; 


<10 


-3 


B-0654, B0738 


64 


D'=' 


YSG-33 


1.02 


<10^'^ 


1.47x10^^ 


0.06 


7.99; 


<10 


-4 


B-014 (V-131) 


65 


B" 




1.94 


<10^^ 


7.66x10^'' 


0.15 


4.11; 


<10 


-4 


C-Z242 (V-486) 


66 


B" 




5.25 


< 10''"' 


2.95x10^^ 


0.21 


1.58; 


<10 


-4 


H-931, H945, H-944, H-939 


67 


C" 


YSG-19 


9.20 


KlO^'' 


1.30x10^"* 


0.47 


6.93; 


<10 


-4 


B-614, V-100 


68 


B" 




2.49 


KlO^*^ 


4.65x10^^ 


2.27 


2.47; 


<10 


-4 


H-907, H-933 


69 


B" 




5.50 


KlO^*^ 


2.82x10^^ 


1.17 


1.50; 


<10 


-4 


H-823, V-502 


70 


B" 




3.08 


<10"^ 


2.52x10^'* 


0.29 


1.35; 


<10 


-3 


H-793, V-227 


71 






1.59 


<10^^ 


8.05x10^^ 


0.19 


4.31; 


<10 


-4 


C-Z040 



Note. — Col. (1): GMC ID. Col. (2): GMC classifications: Type A shows no signature of massive star formation; Type B is associated 
with Hll region(s); Type C is associated with Hll region(s) and less than lOMyr-old stellar group(s); and Type D is associated with Hll 
region(s) and 10— 30 Myr-old stellar group(s). Col. (3): Young stellar groups (YSGs) associated with the GMCs. Col. (4): Luminosity 
of HcK, as measured over the extent of a GMC. Col. (5): Luminosity of 24 / J.m, as measu red over the ex tent of a GMC . Col. (6): The 
Ha attenuation is calculated as Ahq ^ 2.5 log[l + 0.038 L(24 /im)/L(Ha)] jCalzetti eTa l. 2005 Presco tFet al.|2007[ . Col. (7): Star 

1 Detween the extinction-corrected Ha line emissio 



formati on rate (SFR), derive d from L(24/im) and L(Ha) using the relation between the extinction-co rrected Ha line emission and the 
SFR inlCa l zetti et at ||20 07ll. Col. (8): Name of the associ ated Hll regions : B- XX from fBoulesteix erar|jl974[ ; C-XX from[Courtes| 
|et al.| ) |^?7^ H-A A irom [Hodge et al.| l [T999[ l; Y-XX from [Verley et al.| | |2007[ |. The names within parentnesis indicate a counterpart 
of Hit regions in the 24 ^ra catalog. 

^GMCs at the edge of the field-of- views. The values just provide lower limits. 

'^The selected GMCs with more accurate classification. See details in the text in Section |3.4| 
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A. Spatial Comparison of GMCs, Stellar Groups and Hii Regions for Individual GMCs 

Figures 18 - 78 (to be available online only) show the spatial comparison of GMCs, -R3-2/1-O7 YSGs 
and Hii regions for the identified GMCs except those at the edge of the observed field, GMC-18, GMC-28, 



GMC-47, GMC-51, GMC-58 GMC-71, and also four GMCs that have been mentioned in Section 3.5 
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